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Abstract 

We review the role of two-photon exchange (TPE) in electron-hadron scattering, focusing in 
particular on hadronic frameworks suitable for describing the low and moderate region relevant 
to most experimental studies. We discuss the effects of TPE on the extraction of nucleon form fac- 
tors and their role in the resolution of the proton electric to magnetic form factor ratio puzzle. The 
implications of TPE on various other observables, including neutron form factors, electroproduc- 
tion of resonances and pions, and nuclear form factors, are summarized. Measurements seeking to 
directly identify TPE effects, such as through the angular dependence of polarization observables, 
nonlinear e contributions to the cross sections, and via e^j) to e~p cross section ratios, are also 
outlined. In the weak sector, we describe the role of TPE and 7Z interference in parity-violating 
electron scattering, and assess their impact on the extraction of the strange form factors of the 
nucleon and the weak charge of the proton. 



*Dedicated to the memory of John A. Tjon. 
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1 Introduction 

The electromagnetic probe has been a primary experimental tool in the study of hadron physics for 
many decades. Electromagnetic interactions are extremely well understood, and the pointlike nature 
of electrons and muons make them ideal probes of the internal structure of hadrons. Lepton scattering 
experiments provide the cleanest information available on fundamental quantities such as hadron form 
factors and parton distributions. In addition, the relatively small value of the electromagnetic coupling 
means that measurements on nuclei probe the entire nuclear volume, without the significant attenuation 
of the beam over the length of the nucleus that yields a dominance of surface effects in some hadron 
beam measurements. 

Because of the power of the electromagnetic probe, a great deal of our information on the structure 
of the nucleon comes from unpolarized measurements of the inclusive lepton-nucleon cross section. More 
recently, polarized beams have been used to provide additional information on the spin structure of the 
nucleon and to improve our knowledge of the nucleon form factors. As one of the most fundamental 
observables characterizing the composite nature of the nucleon, electromagnetic form factors have over 
the past few decades provided considerable insight into the nucleon internal structure, with the electric 
(Ge) and the magnetic (Gm) form factors encoding the (transverse) spatial distributions of the nucleon's 
charge and magnetization (for reviews see Refs. [H El E])- While these can in principle be extracted 
from the unpolarized cross sections, polarization measurements have played a critical role in studies of 
the nucleon form factors over the last decade. 

The polarization measurements have led to a renaissance in studies of the structure of the proton 
and neutron, providing significantly improved measurements of their form factors over a wide range of 
momentum transfer (which, loosely speaking, reflects the inverse resolution at which the structure is 
probed). However, while these measurements have considerably improved the precision with which the 
form factors could be extracted, they also revealed a significant discrepancy with extractions from the 
unpolarized cross sections in kinematic regions where both techniques provide precise measurements. 
Because these, and essentially all other electron scattering measurements, are analyzed in the one-photon 
exchange or Born approximation, this discrepancy led to a serious reexamination of the possible role 
played by two-photon exchange (TPE) corrections. Early measurements and calculations suggested 
that TPE effects were small, although recent studies have provided convincing evidence that these 
corrections can nonetheless be extremely important in specific observables. 

In the 10 years since the form factor discrepancy was confirmed, a great deal of progress has been 
made in understanding TPE contributions. There have been several approaches used to calculate TPE 
corrections for a variety of reactions and observables, as well as a significant effort aimed at constraining 
these experimentally. At present there are calculations of TPE that are consistent with all existing 
experimental constraints, which can explain the form factor discrepancy, and which allow extraction 
of the proton form factors without yielding a significant theoretical uncertainty in the extraction. In 
addition, a set of experiments is underway which will allow for direct verification of these calculations, 
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and thus provide a final resolution to the issue. With extensive experimental checks of the TPE 
calculations, reliable estimates can be made for other reactions, and it will be possible to identify other 
cases where TPE effects may be large enough that they pose concerns for the interpretation of precision 
experiments. 

It is therefore timely to review the recent experimental and theoretical efforts dedicated to studying 
TPE in electromagnetic processes, as well as the impact of two-boson (7 or Z) exchange corrections 
which enter into weak processes induced by electromagnetic probes. In Sec. |2] we provide an overview 
of the relevant electron scattering formalism, including definitions of kinematics and scattering am- 
plitudes. Section [3] summarizes the initial evidence for TPE effects, as well as constraints from early 
measurements. The formalism and calculations of TPE corrections are introduced in Sec. HI and their 
impact on elastic scattering measurements examined in Sec. [51 Implications of TPE for observables 
in other electron-hadron scattering reactions are discussed in Sec. [61 Finally, in Sec. [7] we review the 
role of two-photon and 7Z interference in parity-violating electron scattering and their impact on the 
extraction of the strange form factors of the nucleon and the weak charge of the proton. We end with 
conclusions and the outlook for future studies of two-photon exchange in Sec. [H 

2 Elastic electron— nucleon scattering 

In this section we define the general kinematics of elastic electron-nucleon scattering (Sec. 12.11) . and 
present amplitudes and cross sections in the one-photon exchange or Born approximation (Sec. 12. 2p . 
Following this we discuss the extraction of the electromagnetic form factors in the Born approximation 
using the Rosenbluth separation and polarization transfer methods (Sec. 12. 3p . 

2.1 Kinematics 

For the elastic scattering process eA^ — eA^ the four-momenta of the initial and final electrons are 
labeled by k and k', with corresponding energies E and E', and of the initial and final nucleons by 
p and p', respectively. The four-momentum transfer from the electron to the nucleon is given by 
q=p'— p = k — k', with = — > 0. Conventionally the scattering cross section is defined in terms 
of and the electron scattering angle 9, or equivalently the dimensionless quantities 

4M2' iy2^^^i^^y 
where v = k ■ p/M"^ — r. In the target rest frame the variable e is related to the scattering angle 9 by 

e = (^l + 2(l + r)tan2 , (2) 

and is identified with the relative flux of longitudinal virtual photons. In terms of r and e the incident 
electron energy is 

E = M (r+v/r(l + r)(l+ £)/(!- 5)) , (3) 

and the scattered electron energy is E' = E — 2Mt. 

One can also express the elastic cross section in terms of any two of the Mandelstam variables s 
(total electron-nucleon invariant mass squared), t, and u, where 

s = (k + pY = {p' + k')\ t={k- k'f = -Q\ u = {p- k'f = {p' - k)\ (4) 

with the constraint s + 1 + u = 2M'^ + 2ml. Furthermore, the variable u is related to the Mandelstam 
variables by z/ = (s — u)/{4M'^). The electron mass rrie can generally be ignored at the kinematics of 
interest here. In particular, there are no mass singularities in the limit mg — )■ in either the one-photon 
exchange amplitude or the total TPE amplitude. 
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Figure 1: Elastic electron-nucleon scattering in the one-photon exchange (Born) approxi- 
mation. Particle momenta are indicated in parentheses. 



2.2 Born approximation 

In the Born approximation (see Fig. [1]), the electron-nucleon scattering invariant amplitude can be 
written as 

= --^J^i^Ji^^ (5) 

where e is the electric charge, and the matrix elements of the electromagnetic leptonic and hadronic 
currents are given in terms of the lepton {ue) and nucleon (mat) spinors by 

i-f^i = Ue{k')-f^Ue{k) , Ji^ = UN{p')T^f{q)uN{p) ■ (6) 

Our metric and other unstated conventions follow Ref. [1]. Note that other conventions for amplitudes 
have also been used in the TPE literature [SI El El El E] • The electromagnetic hadron current operator 
is parametrized by the Dirac [Fi) and Pauli {F2) form factors as 

r^;(g) = rF,{Q') + -i^F^m. (7) 

where M is the nucleon mass. In terms of the amplitude M.^, the differential Born cross section is given 

by 



[aMQ^EJ ~ e{l + T) 



\M^\ = — — -Or. 



where a = c'^/Att is the electromagnetic fine structure constant, and the Mott cross section for the 
scattering from a point particle is 

a^E' cos\e/2) 
4^3 sin^(e/2) 
The reduced Born cross section is given by 

an = eGUQ') +TGliQ'), (10) 

where the Sachs electric and magnetic form factors Ge,m{Q'^) are defined in terms of the Dirac and 
Pauli form factors as 

= F,m-TF,m, = F.m+F^m- (n) 

The form factors are normalized such that G''^J^\o) = 1 (0) and G^}"''(0) = /ip(n) for the proton (neu- 
tron), where fj,p(n) = 2.793 (—1.913) is the proton (neutron) magnetic moment. 
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2.3 Form factors in the Born approximation 

The standard technique to extract the electric and magnetic form factors of the proton has been the 
Rosenbluth, or longitudinal-transverse (LT), separation method [TU]. Using the fact that the Born 
level form factors in Eq. ( ITOj) are functions of only, analyzing the cross section as a function of the 
longitudinal photon polarization e at fixed allows one to extract G\j from the e- intercept, and G"^ 
from the slope in e, once standard radiative corrections have been applied. The cross section at ^ = 180° 
{e — )■ 0) depends only on the magnetic form factor Gm, while the cross section at smaller angles is a 
combination of magnetic and electric contributions. Because of the e/r weighting of G"^ relative to 
G^, the contribution from the electric form factor to the cross section is suppressed at large Q^. The 
proton form factor ratios extracted via the Rosenbluth technique have generally been consistent with 
scaling, \Ge\ ~ IGm/a^pI IHl Ull HSl [H]. Note that because the cross sections are sensitive to 
the squares of the form factors, the signs on the form factors cannot be determined from Rosenbluth 
separations alone. 

An alternative method of extracting the ratio R utilizes polarization degrees of freedom to increase 
the sensitivity to the electric form factor at large Q^. Here, longitudinally polarized electrons are 
scattered from an unpolarized proton target, with the polarization of the recoiling proton detected, 
ep — )■ ep. The polarization of the incident electron (or recoil proton) is characterized by the spin 
four- vector [51 [15] 

= (^-X + k-^^), (12) 



m m{m + E) _ 

where m and E are the particle's mass and energy, and the three-dimensional spin vector specifies the 
spin direction in the rest frame. In the limit — )■ 0, the spin four- vector — )■ (0; <^). Since is a unit 
vector, C,"^ = 1, from Eq. ([12]) one has = —1 and k ■ s = 0. For incident electron energies E ^ me, 
the electron spin four-vector s(f can be related to the electron helicity h = (^^ ■ k hj 

^ h — . (13) 
me 

The coordinate axes are chosen so that the recoil proton momentum p' defines the z axis, in which 
case for longitudinally polarized protons one has C,p = P - the Born approximation the elastic cross 
section for scattering a longitudinally polarized electron with a recoil proton polarized longitudinally is 
then given by 

h CTMott — T7—\ -;-, — tan - Gm . (14) 



dn '""^^ M V 1 

For a proton detected with transverse polarization the x axis is defined to be in the scattering plane, 
X = y X z, where y = k x k' defines the direction perpendicular, or normal, to the scattering plane. 
The cross section for producing a transversely polarized proton, Cp ' P' = 0, is given by 

/i CTMott 24/— - — tan- GeGm, (15) 



dn "'"'^ y i + T 2 

while in the Born approximation the normal polarization is identically zero. Taking the ratio of the 
transverse to longitudinal proton cross sections then yields the ratio of the electric to magnetic proton 
form factors, 

/r(l + £) Pt E + E' 9 Pt Ge 

where Pl and Pt are the polarizations of the recoil proton longitudinal and transverse to the proton 
momentum in the scattering plane, proportional to the longitudinal and transverse cross sections in 
Eqs. ( !T^ and ([T5]l . respectively. 
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Figure 2: Ratio of proton electric to magnetic form factors as extracted using Rosenbluth 
(LT) separation [11] (squares) and polarization transfer measurements |T6l [18] (circles). 
Figure adapted from Ref. 



In a series of recent experiments at Jefferson Lab |l6l[I7l[T8lll9l[2nil2ll|22l[23l[2lll25], the polarization 
transfer (PT) technique has been used to accurately determine the ratio Ge/Gm up to = 8.5 GeV^. 
In addition, there have been complementary measurements using polarized targets at MIT-Bates [26] 
and Jefferson Lab [27]. The results, illustrated in Fig. |2l are in striking contrast to the ratio obtained 
via LT or Rosenbluth separations, showing an approximately linear decrease of R with which is in 
strong violation of the scaling behavior (see also Refs. [H [21 [281 129]). 

The discrepancy between the LT and PT measurements of Ge/Gm has stimulated considerable 
activity, both theoretically and experimentally, over the past decade. Attempts to reconcile the mea- 
surements have mostly focused on improved treatments of radiative corrections, particularly those 
associated with two-photon exchange, which can lead to additional angular (and thus e) dependence 
of the cross section. In the following sections we discuss experimental efforts to better understand 
the discrepancy, and then describe theoretical attempts to compute TPE corrections and assess their 
impact on various observables. 



3 Experimental observables and measurements 

3.1 Verification of the discrepancy 

The striking difference between Rosenbluth [30] and the early polarization transfer [161 [H] measure- 
ments of the proton electromagnetic form factors shown in Fig. [2] led to significant activity aimed at 
understanding and resolving this discrepancy. It was noted early on [TH] that there was significant scat- 
ter between the results of different Rosenbluth extractions [HI [311 [321 [331 131], as illustrated in Fig. [3l 
suggesting that the problem was related to the cross section measurements. At high Q^, Ge yields only 
a small, angle- dependent correction to the cross section, leading to the possibility that a systematic 
difference between small- and large- angle measurements could yield large corrections to Ge/Gm that 
would increase in importance with increasing Q^. Thus, it was initially suggested that the observed 
difference may be due to some experimental error in one or more of the cross section measurements, 
rather than a true discrepancy between the techniques. 

Because most extractions combine multiple data sets covering different angular regions, uncertainty 
in the relative normalization of the data sets leads to large uncertainties in the form factors that can 
be strongly correlated between points at different values. Many extractions allow the normalization 
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Figure 3: Ratio of proton electric to magnetic form factors as extracted from the initial 
high polarization transfer measurement [16] (stars) and previous Rosenbluth (LT) sepa- 
rations [III EH I 



of each experiment to vary as part of the fit, to improve consistency in the overlap regions, but the 
uncertainty in the determination of the relative normalization factors was usually neglected. At high 
Q^, a change in the normalization between large- and small- angle measurements can systematically 
shift all of the high results within the extraction. 

A detailed reanalysis of the world's cross section data for moderate to high values showed 
that the cross sections were consistent when accounting for the quoted normalization uncertainties 
in the different data sets, except for some of the small-angle data from Ref. [TT]. Excluding these 
data and performing a global extraction of G^; and Gm yielded results that were still in significant 
disagreement with the polarization transfer extractions. Furthermore, it was not possible to reproduce 
the polarization results by varying the normalization of the data sets within their quoted uncertainties, 
or by simply excluding one of the 20 data sets included in the analysis. However, it was noted that a 
systematic correction with an e dependence of ~ 5% for Q"^ > 2 GeV^ would be able to explain the 
discrepancy between the polarization measurements and the Rosenbluth extractions [T2] . 

Shortly thereafter, several extremely high precision measurements of Ge/Gm were made using a 
modified Rosenbluth technique [TU [35]. Because the extraction is very sensitive to angle-dependent 
corrections, the measurement used detection of the recoil proton, rather than the scattered electron. In 
the traditional Rosenbluth separation, varying e while keeping fixed means that the scattered electron 
energy and angle vary with e, as does the scattering cross section. This means that any momentum- or 
rate-dependent efficiencies or systematic corrections yield e-dependent effects. With proton detection, 
the momentum of the detected particle does not vary for measurements at a fixed value, and the 
cross section for detection of the recoil proton is only very weakly dependent on angle. These and several 
other corrections and systematic uncertainties are thus strongly suppressed in the extraction of Ge/Gm 
[m |35]. In this measurement, care was taken to ensure that the relative cross section uncertainties 
were minimized, even though this could increase the absolute uncertainties, because only the relative 
uncertainties enter into the ratio Ge/Gm- The experiment yielded measurements for Q^=2.64, 3.2 and 
4.1 GeV^ with precision comparable to the recoil polarization results, as shown in Fig. HI This confirmed 
a clear and significant discrepancy with the polarization measurements, which ruled out the possibility 
that the difference was simply caused by a small systematic error in some of the earlier measurements. 

It is worth noting that the radiative corrections applied to the elastic scattering cross sections have 
evolved somewhat over time. While the early work of Mo and Tsai [STJ [38] provides the standard basis 
for applying radiative corrections, the SLAC measurements in the 1990s [TTl[3l] accounted for additional 
radiative correction terms such as vacuum polarization contributions from higher mass leptons, which 
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Figure 4: Proton electric to magnetic form factor ratio ^IpGe/Gm (left) and {^IpGe/GmY 
(right) extracted from the high-precision Super- Rosenbluth experiment [H] (filled circles), 
compared to polarization extractions [T6| IT8] (triangles) and a global analysis of previous 
cross section measurements [36] (crosses). Note that the slope of the reduced cross section in 
the Rosenbluth measurements is directly sensitive to {fJ^pG e / G m)^ , so the right figure best 
shows the significance of the discrepancy at the cross section level. Figure adapted from 
Ref. 



have also been included in other recent extractions. In the global analysis |36], these corrections were 
also included for older measurements to minimize the difference in the radiative correction procedures 
over the full body of data. However, the additional terms have essentially no dependence on the electron 
scattering angle, and therefore have minimal impact on the extraction oi Ge at large Q^. 

3.2 Early estimates of two-photon exchange corrections 

Once it was clear that there was a systematic discrepancy between the Rosenbluth and polarization 
transfer techniques, it was important to determine how large of an effect would be needed to explain 
the discrepancy, and determine if TPE contributions of this size could be ruled out by existing measure- 
ments. Early attempts to calculate the TPE contributions in the 1950s and 1960s, using only the proton 
in the intermediate state [23110], or including some excited intermediate states [H j W2 \ 1 ^ 1 ^ 1 ^ WTL 
yielded extremely small TPE contributions, typically well below one percent of the Born cross sectioio. 
It was difficult, however, to determine how reliable these estimates were, especially at large values. 
Even before the form factor discrepancy became clear, there were some estimates examining the radia- 
tive corrections [H] and TPE contributions [19] . Coulomb distortion, corresponding to the soft-photon 
contribution to TPE, had also been investigated. This was mainly of interest at low [20], where 
it contributes a significant portion of the full TPE correction. After the observation of the discrep- 
ancy, it was examined at higher [51], but found to have a relatively small effect compared to the 
observed discrepancy. For the most part, investigations have focused on the effect of TPE corrections 
[3 SSI [521 |53l HH |55] beyond the soft-photon approximation. Calculations of TPE corrections will be 
presented in detail in Sec. lU and their impact on the measurements in Sec. |3 

Initial investigations focused on using measurements of the discrepancy to determine the nature of 
the corrections required to explain the data. Analyses of the discrepancy that assume it is due to missing 
strength in the cross section measurements [121 ES] indicated that the difference could be explained by 

^The TPE diagram also contains an infrared- divergent contribution, which is canceled by a corresponding divergent 
term in bremsstrahlung emission, as discussed in more detail in Sec. 14.11 
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Figure 5: The e dependence of the reduced cross section as predicted from the polarization 
transfer resuhs for Ge/Gm (red dashed hne), and as measured by Jefferson Lab experiment 
EOl-OOl [H] (circles; blue solid line is best fit). If the polarization transfer represents the 
true form factors, TPE yields more than half of the e dependence at 2.64 GeV^ (left), and 
85% at 4.1 GeV^ (right). The PT measurements give only the slope of u/j, and the curve 
has been arbitrarily normalized to match the LT separations at e = 0. 



an error in the e dependence of the cross section of approximately 5 — 8% for > 2 — 3 GeV^. Similar 
results are obtained by comparing global analyses of the form factors, with and without the high-Q^ 
Rosenbluth data included [56], [57], under the assumption that the difference will be entirely due to 
the TPE contributions. This approach gives a less direct extraction of the TPE contribution to the 
cross section and is sensitive to other assumptions made in the fit, but it allows the extraction to be 
performed in the context of a more global analysis of form factor data. 

Figure [5] shows the e dependence observed in the Super- Rosenbluth experiment [H] and as expected 
from the polarization transfer measurements [l6l[T8], assuming the Born approximation. The polariza- 
tion transfer measurement constrains only the slope, not the value of the cross section, and has been 
arbitrarily normalized to agree with the Rosenbluth extraction at e = 0. If some systematic correction 
to the cross section explains the difference, then the e dependence shown by the dashed line is related 
to the contribution from the electric form factor, and the remaining slope must come from the missing 
correction. This correction should be close enough to linear such that it does not spoil the linearity 
expected from the Rosenbluth formula, since the reduced cross section is consistent with a linear e 
dependence within the current uncertainties. 

If the discrepancy is entirely related to TPE contributions, then there are additional constraints 
that can be included in such extractions. For example, TPE corrections must be zero at e = 1, as 
discussed in Sec. HJ so one would expect the two lines in Fig. [5] to meet at e = 1 rather than e = 0. Such 
constraints were included in later attempts to extract the TPE contribution to the cross section. Two 
such extractions jSHl EH] use symmetry arguments to choose functional forms for the TPE corrections 
and then fit the difference between Rosenbluth and polarization data to extract the TPE contributions. 
While the symmetry arguments do not give a unique e dependence, they suggest that it may be more 
natural to take the correction to be a simple function of x = \J {1 + e) / {1 — e) rather than a function 
of 9 or e. Note, however, that the specific form chosen in Ref. pS] is divergent for ^ — )■ 0°, and so does 
not yield the correct limit for e = 1. Similar extractions are performed in Refs. [511 EDI ED E2], again 
fitting only the TPE contributions to the cross section, but making somewhat different assumptions 
about the functional form. A detailed comparison of the assumptions and procedures of many of these 
extractions can be found in Ref. [6T| . 
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It is worth noting that many of these TPE extractions are only rehable at high values. All 
of the extractions discussed thus far assume that the TPE impact on polarization measurements is 
negligible. While this assumption is reasonable at high values, where the Rosenbluth extraction of 
is extremely sensitive to TPE, it is not as good at lower values where both the polarization 
and cross section measurements may have small but similar TPE contributions. This can introduce 
additional uncertainty in the form factors at low in cases where both the TPE contribution and 
form factors are extracted. A clear example is Ref . [36] , where the TPE contribution extracted at large 

is applied as a correction at all values in extracting the form factors, leading to an overestimate 
of TPE effects at low Q^. 

The first attempt to use the discrepancy to extract the TPE amplitudes was made in Ref. [52], 
which took into account the impact of TPE on the polarization measurements as well as on the cross 
section data. Because there are three TPE amplitudes, each of which depends on and £, significant 
assumptions have to be made about their e dependence. Thus, different extractions [52], |63l |63j yielded 
very different amplitudes, depending on the constraints or approximations used to select the amplitudes 
and functional forms used in the analysis. However, because the dominant effect at high is related to 
the TPE contribution to the cross section, all of these analyses yield similar results in this region: an e 
dependence of 5-8% at high Q^, although with different overall normalizations. Later analyses [HS] ES] 
were able to make less model-dependent extractions with the inclusion of new information on TPE in 
polarization measurements [67] . These results will be discussed in more detail in Sec. 15.31 

Early calculations were performed in hadronic [5] and partonic [55] frameworks, yielding corrections 
which could explain roughly half of the discrepancy at large values. Calculations at the parton 
level in the double logarithm approximation [51] yielded a different form for the e dependence, with 
nonlinearities appearing at large e. Finally, invariance under C-parity and crossing symmetry were 
used [53] to argue that the TPE corrections should depend on x = + £)/(! — e). The range of 
predictions of these models is shown in Fig. [6] for ~ 2 — 3 GeV^, where each model has been scaled 
to give an overall e dependence of the size needed to resolve the discrepancy between the Rosenbluth 
and polarization measurements. For the partonic model [55] , the calculation is not expected to be valid 
at low or e values. Many of these calculations have been updated since, and the most recent results 
will be discussed in Sec. IH but these were the models available at the time that measurements to further 
examine the TPE corrections were being considered. 
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Figure 6: The TPE correction 5, defined via a = ao{l + S), where (Tq is the Born cross 
section, to the elastic ep cross section from the calculations of Rekalo et al. [53] (dashed), 
Chen et al. [55] (dot-dashed), Blunden et al. [5] (solid), and Afanasev et al. [51] (dotted) for 
~ 5 GeV^, after scaling the calculations to yield an e dependence of approximately 6% 
over the e range of existing data. 
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While scaled to yield approximately the same e dependence, the different models would nonetheless 
yield significantly different results for the extracted form factors. The Chen et al. calculation ^55j has 
little e dependence for e > 0.5, meaning that it would have little impact on Rosenbluth extractions if 
they did not have data at low e. The form of the correction in Ref. [53] yields large effects for e — )■ 1, but 
small corrections to G]\i for measurements at low e. The other calculations have the largest corrections at 
low e, and thus larger corrections to direct measurements of Gm- In addition, some calculations predict 
large deviations from the linear e dependence of the Born approximation at very high e, while others 
predict the greatest deviations at low e. Based on this large range of predictions, several experimental 
tests, using existing data or new measurements, were performed to provide independent constraints on 
any possible TPE contributions. 

3.3 Experimental signatures of two-photon exchange contributions 

Two-photon exchange contributions to elastic electron-proton scattering manifest themselves in several 
different ways. The real part of the TPE amplitude modifies both the unpolarized cross section and 
the polarization transfer components used to extract Ge/Gm, and thus measurements sensitive to the 
real part of TPE amplitudes yield constraints that are directly relevant to the form factor discrepancy. 
These contributions, however, must be disentangled from the dominant single-photon contributions. The 
imaginary part of the amplitudes, on the other hand, leads to non-zero values for the Born-forbidden 
normal asymmetries (Sec. 15. 5p . which allow two-photon effects to be directly isolated. However, since 
these are not directly connected to the form factor discrepancy, they were not the initial focus of 
experimental investigations. We summarize here some of the early experiments that provided constraints 
on the TPE contributions. 

3.3.1 Comparisons of positron— proton and electron— proton scattering 

In unpolarized electron-proton scattering the TPE contributions cannot in practice be separated from 
the Born cross section empirically. The cleanest means of identifying TPE effects in unpolarized ep 
scattering is to compare positron-proton and electron-proton cross sections, where the interference 
between one- and two-photon exchange has the opposite sign for positron and electron beams (Sec. 14. 1| ). 
While experimentally difficult, these were considered crucial tests of the electron scattering technique, 
and several experiments were performed in the 1960s to study possible TPE effects. The comparisons of 
e+p and e~p |6H1 EH [TOl [TH [721 [731 El] and fi~^p and fi^p [75j scattering were interpreted as supporting 
the conclusions of the early calculations that two-photon corrections were extremely small (< 1%). 
Figure [7] (left panel) shows a compilation of all such measurements for elastic scattering as a function 
of Q^. However, the low intensity of e"*" (/i"^) beams made precise measurements nearly impossible for 
large or small e values. 

The data were reexamined ^76j in light of the form factor discrepancy, which suggested an e- 
dependent effect. Measurements at high were limited to small scattering angle (large e) and so 
could not set a meaningful limit on angle-dependent TPE contributions if they were small at forward 
scattering angles. Figure [7] (right panel) shows the ratio of positron to electron cross section as a func- 
tion of e for measurements at < 2 GeV^. There is some evidence of a charge- dependent term in the 
e^p elastic cross section ratio at small values of e, although the data at low e are not very precise. A 
linear fit yields a 6% e dependence in the positron to electron ratio, implying a 3% e dependence for the 
electron-proton cross section. This is roughly half of the e dependence needed to explain the discrep- 
ancy at high Q^, but the average value of the points showing nonzero TPE is less than 0.5 GeV^, 
well below the region of the observed discrepancy. In addition, these measurements need to apply 
a model-dependent correction for the charge-dependent bremsstrahlung contributions to isolate TPE 
effects, and such corrections have not been applied in a consistent fashion for all of the measurements. 
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Figure 7: The ratio of positron-proton to electron-proton scattering as a function of 
(left) and e (right). The e dependence plot includes only measurements below = 2 GeV^. 
The solid curve is a linear fit to the data yielding a slope of —0.057(18). Figure adapted 
from Ref . 



Because of the limited statistics, the low values of the large angle measurements, and the model 
dependence in extracting the TPE contributions, these data are insufficient to make a strong conclusion 
about the presence of TPE effects. There are three experiments which aim to improve the precision and 
kinematic coverage of e^p comparisons [771 [THl [79] . Two of the measurements [771 [78] have completed 
data taking, and the third is expected to begin in the near future. While all of these will be limited to 

<i 2 GeV^, they will make precise measurements over a significant e range, and allow direct tests 
of calculations of the TPE effects. Additional details on the planned measurements are presented in 
SeclO 



3.3.2 Improved measurements of the Rosenbluth— polarization transfer discrepancy 

Given the difficulty of making precise comparisons of positron and electron scattering at high and 
large scattering angle, one must complement these direct measurements with other studies that can be 
used to constrain TPE contributions. The most compelling evidence to date is the discrepancy between 
Rosenbluth and polarization measurements of the proton form factor ratio Ge/Gm which can be used to 
constrain the overall size of the e-dependent TPE effects under the assumption that TPE fully explains 
the discrepancy and that all other corrections are accounted for. 

Even at high values, the discrepancy is only at the 2-3 a level when examined as a function 
of {^IpGe/GmY, Fig. [D (right panel), which is directly related to the measured slope of the reduced 
cross section. While new or updated high polarization measurements are now available [221 121], the 
comparison is clearly limited by the uncertainties in the Rosenbluth separation measurements. Jefferson 
Lab experiment E05-017 [80] made an extended set of Super- Rosenbluth measurements, using the same 
technique as Ref. [H], but covering a much larger range, and the data are currently being analyzed. 

Below = 1 GeV^ the sensitivity to TPE corrections is smaller, and significant improvements 
in both the polarization and Rosenbluth measurements are necessary to search for indication of TPE 
contributions. Extended Rosenbluth measurements from Mainz [81] and high-precision polarization 
measurements from Jefferson Lab [231 123 E2] are now available. The new data show reasonable agree- 
ment between Rosenbluth and polarization measurements below = 1 GeV^, suggesting that TPE 
contributions may be small in this region. However, because the impact of TPE corrections is reduced 
at low Q^, and the Rosenbluth extractions [81] did include a partial implementation of TPE (using the 
Coulomb distortion correction in the = limit from Ref. [39j), it is not yet clear to what extent the 
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agreement sets significant limits on TPE. 



3.3.3 Experimental limits on nonlinearities 



In the Born approximation, the reduced cross section at fixed depends linearly on e. Any deviation 
from linearity must come from terms that are not included in the standard radiative correction proce- 
dures. Thus, the difference between polarization and Rosenbluth measurements of the form factors is 
related to the average linear contribution of TPE, while any deviation from linearity is a clear indicator 
of effects beyond the Born approximation. At lower values, observing such a deviation would be a 
clear signature of effects beyond one-photon exchange and would provide quantitative information on 
the nonlinear component of the such effects. At large Q^, the contribution from Ge becomes small 
enough that almost all of the e dependence comes from TPE (or other corrections), as shown in Fig. |5l 
In this region, the e dependence of the reduced cross section allows one to isolate the full e dependence 
of the contributions beyond the Born approximation. 

A detailed examination of nonlinear contributions was performed by fitting the reduced cross section 
to the form 

(jR = Po[l + Piis - 0.5) + P2is - 0.5)2] _ (^7) 

Note that the quadratic term is expanded around e = 0.5 such that P2 represents the relative magnitude 
of the nonlinear term compared to the average cross section, rather than the e = cross section, which 
becomes extremely small at low Q^. The coefficient P2 provides a simple measure of the relative size 
of nonlinear terms, and the uncertainty 6P2 can be used to set limits on 0{e'^) terms. Conventional 
Rosenbluth separation measurements have found P2 to be consistent with zero, and the best constraint 
[3l] yields 6P2 ~ 10.5%. The recent Jefferson Lab EOl-OOl data [l4j obtained improved limits on P2 
{SP2 = 4.4%) by detecting the struck proton rather than the scattered electron. 
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Figure 8: The e dependence of the reduced cross section from SLAC NEll [34J (using 
only data from the 8 GeV spectrometer) (left) and for Jefferson Lab experiment EOl-OOl 
[13] (right). The solid line is a linear fit, while the dotted lines are quadratic fits with 
P2 = ±0.105 for NEll, and P2 = -0.008 ± 0.044 for EOl-OOl {i.e. la variations on the 
central value). 



Figure [8] shows the reduced cross section as a function of e for the best conventional Rosenbluth 
separation (SLAC NEll [33]), and from the Super- Rosenbluth extraction in Jefferson Lab experiment 
EOl-OOl [13]. The data are consistent with no nonlinear contributions. A global analysis of nonlineari- 
ties, including both elastic scattering, resonance region measurements, and deep-inelastic scattering was 
performed in Ref. [83j. This analysis concluded that all of the measurements are consistent with a linear 
e dependence. For the elastic measurements, the limits on nonlinear contributions, shown in the left 
panel of Fig. [9l are not very significant except for = 2 — 4 GeV^, where there are recent high precision 
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Figure 9: (Left) Nonlinearity parameter P2 (see Eq. ([8])) from a global analysis |83j of 
elastic cross section data. The dotted line indicates the global average, (P2) = 0.019 ±0.027. 
(Right) Weighted average of the deviation of the cross section measurements from linear 
fits, = ((Tdata — crfit)/o"fit- The curve is a quadratic fit to the residual deviation yielding 
a quadratic term of 0.9 ± 2.0%. Figure adapted from Ref. [83] . 



measurements [13] • The right panel shows the difference between the cross section measurements and 
the linear fits, averaged over all Q^, which gives an indication of how well the possible deviations from 
linearity are constrained as a function of e. The improved measurements of Jefferson Lab experiment 
E05-017 should allow for improved constraints over a range in Q^. 

Similar analyses were performed by other groups for ep [58] and ed elastic scattering [84J . In Ref. [58] , 
the nonlinear contributions were parametrized using constraints from C-parity and crossing symmetry. 
The parametrization of the TPE contributions there (relative to the Born cross section) was taken to be 
proportional to x = ^J{1 + e) / {1 — e), which diverges for e — )■ 1. Thus, any data with strong constraints 
at high e yields a very tight constraint on the coefficient of the divergent TPE contributions, but still 
yields large nonlinear contributions at higher e values. While this parametrization is consistent with 
the symmetries being considered by the authors, it is inconsistent with TPE contributions, as discussed 
in Ref. [59], and does not provide meaningful constraints on nonlinear TPE effects. 

Based on similar considerations, the analysis of TPE effects in ed scattering by Rekalo et al. [84] 
expands the TPE contributions as a Taylor series in x. However, since x is always larger than unity, 
and in fact diverges as £ — )■ 1 (^ — )■ 0), the expansion of the TPE contribution in powers of x does not 
provide a meaningful way to examine the data. 

In Ref. [85] estimates were provided for P2 based on the partonic TPE calculations [551 ES] (see 
Sec. 14. 6p . This work found that the calculated nonlinearities were of approximately the same size as the 
limits set by the global analysis [S3], but had a large dependence on the model used for the calculation. 

3.4 Conclusions from early experiments 

While examinations of previous data and the early measurements aimed at understanding the form 
factor discrepancy did not provide clear answers, they were consistent with the idea that larger than 
expected TPE contributions could explain the observed effects. The effects had to be small at forward 
scattering angles, to be consistent with comparisons of positron and electron scattering, and yield an 
angle dependence of ~ 5% at high Q^, approximately linear in e, to resolve the discrepancy. While 
these studies did not provide definitive evidence that TPE were responsible, they still provide some of 
the most significant constraints on certain aspects of the TPE contributions. 

Because these corrections were larger than previous estimates of TPE, significant effort was put into 
understanding them theoretically and placing experimental constraints through new measurements. 
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These efforts were not just limited to examining TPE corrections to the unpolarized cross sections; if 
TPE yields effects on the cross section on the order of a few percent, it may also enter into other ob- 
servables at the same level and contribute to reactions beyond elastic ep scattering. While experimental 
tests of TPE are for the most part limited to elastic ep scattering, measurements of several observables 
can be performed. With a range of measurements, one hopes to constrain models well enough that they 
can provide reliable, semi-quantitative estimates for TPE contributions in other reactions. This will 
allow one to determine where residual TPE effects may have a significant impact in other reactions. 

As mentioned in Sec. 13. 3[ the initial experimental investigations focused on observables directly 
related to the cross section measurements. The initial work focused on existing data, while new mea- 
surements were being performed to improve the precision and kinematical coverage of Rosenbluth sep- 
arations, polarization measurements, and comparisons of positron and electron scattering. Additional 
measurements were proposed to study the e dependence of the polarization observables |67j (Sec. 15. 3p . 
which are e independent in the Born approximation, as well as observables which are forbidden in 
the Born approximation (Sec. 15. 5p . Before addressing these additional measurements, in the following 
sections we discuss calculations of TPE corrections within several theoretical frameworks (Sec. H]), com- 
parisons of these calculations to additional measurements (Sec. [5]), and estimates of TPE contributions 
in reactions beyond elastic ep scattering (Sec. [6]). 

4 Two-photon exchange 

In view of the failure to understand the discrepancy between the Rosenbluth and polarization transfer 
measurements of Ge/Gm in terms of standard radiative corrections, the focus soon turned to revisiting 
the methodologies used in computing the box and crossed-box two-photon exchange corrections, illus- 
trated in Fig. [TOl An experimental-theoretical working group was established at Jefferson Lab, with 
the goal of identifying possible directions for resolving the discrepancy. The first quantitative calcu- 
lation resulting from this renewed focus was made by Blunden et al. 0, who computed the effect on 
Ge/Gm from TPE, incorporating explicitly the nucleon's substructure. Thereafter followed a number 
of other studies, examining TPE in a variety of frameworks, and exploring reactions beyond elastic ep 
scattering. In a parallel effort, Guichon and Vanderhaeghen [52] provided a generalized formalism for 
elastic scattering, allowing for possible TPE contributions, and demonstrated that it was natural to 
have TPE contributions which could significantly change the LT extraction oIGe with minimal impact 
on the PT measurements. 

In this section we review these efforts, paying particular attention to the conventional hadronic-level 
calculations which are most applicable to data analysis at low to moderate values. Before turning 
to the (model-dependent) TPE contributions, we first discuss some general properties of radiative 
corrections in order to set the stage for the recent improvements. 




Figure 10: Box and crossed-box two- photon exchange contributions to elastic electron- 
nucleon scattering. The overall four- momentum transfer to the nucleon is g = k — k' = qi+q2- 
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Figure 11: The complete set of diagrams to order a^, including the virtual (left) and 
bremsstrahlung (or inelastic) (right) contributions. Figure taken from Ref. [l9j . 



4.1 General properties of radiative corrections 

At order a^, radiative corrections to the one-photon exchange cross section cr^, Eq. ([8]), include processes 
arising from the exchange of a second virtual photon ((5virt)) and inelastic scattering processes involving 
the emission of a real bremsstrahlung photon ((5brcm), so that 

CTR = 0-^(1 + Svirt + Kem) ■ (18) 

In analyzing the virtual corrections for a ep scattering, it is convenient to separate terms into "soft" 
parts, which are independent of hadronic structure, and "hard" parts, which are model dependent. Soft 
here implies that any interaction of the second virtual photon with the proton occurs with vanishingly 
small momentum transfer. The soft parts are therefore the same for protons as they are for scattering 
from pointlike particles [e.g. in e~fi'^ scattering). All of the infrared (IR) divergences for the virtual 
diagrams are contained in the soft parts, and cancel in the total amplitude. 

If we denote the amplitude for all one-loop virtual corrections, illustrated in Fig.[TTl by A^i_ioop, then 
M i-ioop can be written as the sum of a "factorizable" soft term, proportional to the Born amplitude 
A^^, and a non-factorizable hard part A^hard, 

-Ml-loop = f{Q^,e)My + Alhard, (19) 

with f{Q^,s) a purely kinematic factor. Therefore ^^irt is given by 

27^e{^^;^^hard} 



\M 
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All of the virtual processes in Fig. [TT] contribute to the soft terms in 5soft- In practice these terms 
dominate, and are the ones accounted for in the standard radiative corrections of Mo and Tsai [2Z1 ES]- 
Furthermore, the functions f{Q'^,e) for the vacuum polarization, self-energy, and vertex corrections are 
e- independent (although the vertex terms are IR-divergent), and therefore have no relevance for the 
LT separation aside from an overall normalization factor. Hence, of the factorizable terms, only the 
TPE can contribute to the e-dependence of the cross section (Fig. [10]). It should be noted that the 
decomposition of the TPE effect into soft and hard parts is not unique, as discussed below. 

The terms which depend on hadronic structure are contained in A^hard, and arise from the proton 
vertex and TPE corrections. For the proton, the hadronic vertex correction was analyzed by Maximon 
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and Tjon |39], and found to be < 0.5% in magnitude for < 6 GeV^. Since the proton vertex 
correction does not have a strong e-dependence, it will not affect the LT analysis, and can be safely 
neglected in examining the form factor discrepancy. 

For the inelastic bremsstrahlung contribution, illustrated in Fig. [TTl the amplitude for real photon 
emission can also be written in the form of Eq. ( I19p . In the soft photon approximation (SPA) one again 
keeps only the factorizable terms that are independent of hadron structure, in which case the cross 
section is a simple kinematic factor times the Born cross section, Eq. ([8]). A significant e-dependence 
arises from this contribution due to the frame dependence of the angular distribution of the emitted 
photon. 

Typically an infinitesimal photon mass A is introduced in the photon propagator as a bookkeeping 
device to regulate the IR divergences. For the box and crossed-box TPE diagrams in Fig. [TOl this 
A dependence is canceled by the bremsstrahlung interference contribution with a soft-photon emitted 
from the electron and proton {i.e. by cutting one of the (soft) photon propagators). This produces 
a correction to the cross section depending on the proton charge Z. A similar cancellation occurs 
for the IR divergences in the electron vertex and proton vertex corrections, which depend on and 
Z'^, respectively. For positron-proton scattering one can take Z = —1, so that only the TPE (and 
bremsstrahlung interference) terms are relevant in comparison with the electron-proton case. 

Although the standard radiative corrections are model independent, there have been significant 
improvements to the original work of Mo and Tsai [STJ [38] , removing many of the mathematical and 
other approximations used in that work. An important contribution was the calculation of Maximon 
and Tjon [IH], whose analysis differs from that of Refs. EB] in two substantive aspects. First, 
they evaluated the inelastic bremsstrahlung cross section without any approximation, using integrals 
given in closed form by 't Hooft and Veltman |87]. The exact expressions are simpler in form than 
the approximate ones given in Ref. [38]. In the limit M — > oo, corresponding to a static Coulomb 
potential, they reproduce exactly the result first given by Schwinger |SB]. Second, in the evaluation of 
the contribution of the TPE diagrams, they make a less drastic approximation than that made in [27] . 
Specifically, in the integrands corresponding to the relevant IR-divergent amplitudes, they make a soft- 
photon approximation for the matrix elements of the current appearing in the numerator (as in [37]). 
but not for the propagators, which appear in the denominator. The required integrals (scalar four-point 
functions) have also been given in [HZIISS]; the resulting expressions are again considerably simpler than 
those obtained in [37], where the soft-photon approximation is also made in the denominators. 

Other improvements have also been made to the work of Mo and Tsai [M|- These include higher 
mass vacuum polarizations terms (not just e~^e~) from Ref. [H], improved implementation of multi- 
photon exchange and angle-dependent bremsstrahlung for coincidence reactions [l8l [90] , and improved 
treatment of multi-photon exchange in the structure function method [91] . The first three were simply 
improvements to the Mo- Tsai framework, made independently of the form factor discrepancy. The last 
[?T] was an attempt to explain the discrepancy. While they found differences between their results and 
Mo-Tsai based radiative corrections, the comparison was not made for identical conditions, and more 
work is required to determine whether there is really any substantial difference. 



4.2 Soft two-photon exchange effects 



The improvements in the computation of the radiative corrections made by Maximon and Tjon 
have a number of important effects on the corrections of order Z. The form of the bremsstrahlung 
correction depends on whether the scattered electron or the recoil proton is detected. For the detection 
of the scattered electron, the bremsstrahlung correction in the SPA is found to be 
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with rj = E/E' the ratio of incident to final electron energies, Li2 is the dilogarithm function, and 
y = 1 + 2r + 2^yT{l + r). Here AE is the maximum allowable energy loss in the lab frame due to 
detector acceptance or experimental cuts, below which one cannot determine that a soft photon has 
been emitted. This exact expression is simpler in form than the approximate one of Mo and Tsai 
[371 EH] • Their expression is too long to reproduce here, but it takes the form 

(^bremlMoT) = InT] \n ( ] + finite terms. (22) 



From this it is clear that the logarithmic dependence on both AE and A in Eqs. fl2T|) and fl22l) is the 
same, and so the difference between the two treatments is a finite kinematic factor. 

The IR-divergent part of the box amplitude can be separated from the IR-finite part by analyzing 
the structure of the photon propagators in the integrand. The two poles, where the photons are soft, 
occur at gi = {q2 = q) and at gi = g (g2 = 0). Evaluating the numerator of the integral at either 
value gives a contribution to the TPE amplitude for the box diagram as 



{s- M^)q^M^ Dois;me,X,M,X), (23) 



2n 

where s = M'^+2ME is the Mandelstam variable in the laboratory frame. The integral over the product 
of four propagators is expressed in terms of the four-point Passarino-Veltman function Dq{s) ^89j, which 
can be evaluated numerically using the program Loop Tools p2] • General expressions for the asymptotic 
expansion of all four-point IR-divergent integrals have been given by Beenakker and Denner [93], who 
find 

D„is; A, M, A) = In (— j In (^ j , (24) 

in the limit (s — M^) ^ ml,meM. Here one is only interested in the real part of this expression. 
The convention for the logarithm is ln{—z) = ln{z) — in when 2; < 0, which is achieved by setting 
z z + iO'^. For the crossed-box one can make use of crossing symmetry. This requires that the 
crossed-box amplitude J^^^°^ obey the relation 

^^^°^(n,t) = +-M5°^(3,t)|^^^. (25) 

where t = q^. Thus the total box plus crossed-box amplitude must be even under the interchange s -H- m, 
or equivalently E f-)- —E'. Note that since — m > 0, the crossed-box amplitude has no imaginary 
part. Combining these expressions, and taking the real part only, the final result for the soft virtual 
correction is [19] 

MMTj) = -^ln7^1n^. (26) 

Note that the dependence on the electron mass has dropped out in the final expression, and the 
logarithmic IR singularity in A is exactly canceled when added to Eq. (|2T1) . 

By contrast, in earlier treatments [33, HSl [9l| the SPA is also applied to one of the propagators, 
for example replacing l/(g — gi)^ by 1/g^ when gi — )■ 0, and 1/g^ by 1/g^ when gi — )■ g. Hence the 
IR-divergent contribution is 

TWlj- = -2^{s-M^)M,^ I Ai ^ 
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Figure 12: Difference between the model-independent corrections of Maximon-Tjon [49] and 
Mo-Tsai [371 ES] for virtual TPE (dashed curves) and virtual -|- real (solid). The upper (lower) 
dashed and solid curves correspond to = 1 (6) GeV^. 



where C^^i^s) is the three-point Passarino-Veltman function [SH], and a factor of 2 accounts for the 
contribution from both poles. To facilitate comparison, note that Tsai [37j introduces the function 
K{—k,p) in evaluating this integral, which is equal to (s — M^)Co(s), with s = {k + p)"^- In this 
approximation one therefore expects 
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for which the crossing symmetry property is manifest. In writing the last line we have again used 
the asymptotic expansion of three-point IR-divergent functions given in Beenakker and Denner [93] . 
Equation ( l28l) is equivalent to the form given by Meister and Yennie [9l], which was used in Refs. |55ll86] . 
Equation (l28ll is larger than Eq. (J26i) by a finite amount of approximately avr - a point that was also 
noted in Refs. [55l |86] . 

Tsai |37] makes a further mathematical approximation, to replace p —p in the box contribution 
only, or equivalently s — )■ — 2ME. Although this additional approximation spoils the crossing 
symmetry, it eliminates the term (^tt)^ that arises in Eq. (1281) from the product of logarithms with 
negative argument, and is therefore closer to the result of Maximon and Tjon, Eq. fl2B]) . At the time. 
Mo and Tsai argued that the resulting expression is closer to the exact calculations of ee scattering 
The final result for the Mo-Tsai soft virtual correction can then be expressed as 
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which is clearly no longer asymmetric under E —E'. 

Because the Mo-Tsai result is the one generally used in existing experimental analyses calculations of 
the full TPE contribution presented here will be taken with respect to 5ir(MoT) rather than (5iR(MTj) 
(except where indicated). It is useful to compare the e-dependence of these two treatments of the soft 
virtual corrections. The difference 5iR(MTj) — (5ir(MoT) is independent of A, and is shown in Fig. [T2l 
(dashed curves) as a function of e for = 1 and 6 GeV^. Note that the difference vanishes as ^ 0. 
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The different treatments of tlie IR-divergent terms already fiave significance for tlie LT separation, 
resulting in roughly a 1% change in the cross section over the range of e. Also shown in Fig. [T2] 
is a comparison for the total correction linear in Z (solid curves), which includes the improvements 
to the bremsstrahlung correction, viz. [5iR(MTj) + 5brem(MTj)] — [5iFi,(MoT) + (5brem(MoT)]. This is 
appropriate when the scattered electron is detected, and is the procedure adopted in the generahzed 
parton distribution calculations of Refs. [551 ES] . From these curves it is clear that most of the difference 
arises from the treatment of the TPE diagrams rather than the bremsstrahlung correction. 



4.3 Hadron structure effects: elastic intermediate states 

Going beyond the Mo-Tsai [371 EH] and Maximon-Tjon [?9] approximations, Blunden et al. [5] pro- 
ceeded to evaluate the box and crossed-box TPE diagrams without resorting to any of the soft-photon 
approximations discussed above. In particular, they considered explicitly the effects of incorporating 
the hadronic structure of the nucleon, parametrized through hadronic electromagnetic form factors. 
Specifically, Blunden et al. computed [5l |6l |95] the total TPE amplitude 



where 



M^^ = A^5°^ + Mff"", (30) 

_^box{xbox) ^ -,e^|^Ljr^^'°^^i^rAF(9i,A) A^(g2,A), (31) 

with the box and crossed-box leptonic tensors given by 

L)ff = u,{k')-f,SF{k-q,,me)-f.u,{k), Lf,°^ = u,{k') Spik - q,,me) J,u,{k). (32) 

Alternatively, one can obtain the crossed-box term directly from the box term by applying the crossing 
symmetry relation Eq. fl25|) . The electromagnetic nucleon elastic hadronic tensor H'^ is given by 

= Ur,{p')T^^{q2)Sp{p + quM)r;{q^)uM. (33) 

Here the electromagnetic current operator Ti^ is given in Eq. (^^, and the fermion (electron and nucleon) 
and gauge boson (photon) propagators are given by 

SF{k,m)= J^^^l. , Ai.(fc, A) = I (34) 

respectively. The infinitesimal photon mass A is introduced to regulate the infrared divergences. One 
can verify explicitly that the integrals in Eq. f l3T|) satisfy the crossing symmetry constraint in Eq. fl25|) . 

The relative correction to the elastic Born cross section, Eq. ([8]), due to the interference of the one- 
and two-photon exchange amplitudes (Figs. [1] and [TOl) is given by 

25e|^^ ^^^^ 

\M.f\ 

Typically, experimental analyses of form factor data apply radiative corrections based on the Mo-Tsai 
prescription [371 EH] or modified approaches based on the same general framework [ll][48j, which include 
approximating the TPE contribution by the IR prescription (5ir(MoT) in Eq. ( I29l) . To determine the 
effect of the full, hadron-structure dependent correction on the data, one must therefore compare 
the e-dependence of the full calculation with that of (5ir(MoT). A meaningful comparison can be made 
by considering the difference 

J = 5^^-5ir(MoT), (36) 
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Figure 13: Difference 5 between the full TPE correction to the elastic cross section [5] and 
the commonly used Mo-Tsai approximation (129]) [23 EH], for = 0.001-1 GeV^ (left), and 
for g2 = i_6 GeV^ (right). 



in which the IR divergences cancel, and which is independent of A. 

The results for the difference S between the full calculation and the MT approximation are shown 
in Fig. [13] as a function of e for several values of from 0.001 to 1 GeV^ (left panel) and 1 to 6 GeV^ 
(right panel). The hadron structure dependent corrections are most significant at low e, where they 
range from 2% to ~ 6% over this range. At the lower values, 5 is approximately linear in e, but 
significant deviations from linearity are observed with increasing Q^, especially at small e. 

Note that assuming that the current operator in Ai^^ has a similar structure off-shell as on- 
shell, with phenomenological form factors at the jNN vertices, introduces model dependence into the 
calculation, as the radiative corrections are used to determine the experimental form factors in the first 
place. In principle this dependence could be removed by iteratively extracting the form factors with the 
inclusion of TPE corrections and feeding those form factors into subsequent TPE calculations. Even 
without iterating, however, because 6^^ is a ratio, the model dependence mostly cancels, provided the 
same phenomenological form factors are used for both Ai^ and Ai-y-y. This was demonstrated in Ref. [6] 
by comparing the results with those obtained using a dipole form 

GniQ') = + ' (37) 

with mass A^, = 0.84 GeV. The model dependence is very weak at the lower value, with virtually 
no effect on the e slope. At larger values the differences increase, but the general trend of the 
correction remains unchanged, so that one can conclude that the model dependence of the calculation 
is quite modest. At ^ 6 GeV^, which will be accessible in future experiments at Jefferson Lab, 
significant deviations from linearity are expected over the entire e range. 

The TPE correction S is also found to be insensitive to the high-Q^ behavior of the Ge/Gm ratio. 
Using form factor inputs from parametrizations obtained by fitting only LT-separated data [521 ES] and 
those in which G^; is constrained by the polarization transfer data [Ml 1^ . the differences are almost 
indistinguishable up to = 6 GeV^ [6j. (Note that the Gm form factor itself also differs by a few 
percent between the various parametrizations.) 

4.4 Inelastic contributions 

While the TPE diagrams in Fig. [10] contain contributions from intermediate states containing all possible 
excitations of the nucleon, the above discussion has thus far been restricted to the nucleon elastic 
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component. In view of the prominent role of the A(1232) resonance, for example, in many hadronic 
reactions, it is important to evaluate its contribution to the TPE amplitude, as well as the possible role 
of other inelastic intermediate states. 

4.4.1 A intermediate states 

The first estimate of the possible role played by excitations of the nucleon in elastic ep scattering was 
by Drell and Fubini [12] . Using dispersion theory methods within a nonrelativistic approximation, they 
found that for energies up to -E ~ 1 GeV, the A(1232) resonance contribution increased the Born 
cross section by ~ 1%. Later, Campbell [HI HS] used relativistic Rarita-Schwinger spinors for the 
A, but retained the soft-photon approximation [38], together with the assumption of Ml transition 
dominance. For large scattering angles, the correction from the A was found to be non-negligible for 

~ 1 — 5 GeV^, but the estimate was limited by the uncertainty in the 'yNA form factor. Partial 
cancellation between the nucleon and A contributions for center of mass energies above 500 MeV was 
observed by Greenhut [16], who used the second Born approximation, neglecting recoil, and assuming 
magnetic only coupling. The conclusion of this work was that the combined and A channels were not 
expected to exceed ~ 1% of the cross section, although the approximations inherent in the calculation 
rendered it less reliable for energies above 1 GeV |46j . 

More recently, Kondratyuk et al. [7] computed the A(1232) resonance contribution within the same 
relativistic framework that was used to compute the nucleon intermediate states in Sec. 14. 3[ without 
invoking any of the approximations of the earlier studies [^ |l5l |16]. Following Refs. [981 EH]! the 
jNA vertex is defined as 

g'^'M + (g^A - ■ pa)] }75, (38) 

where pa and q are the momenta of the outgoing A and incoming photon, with corresponding Lorentz 
indices a and /z, respectively. The factor a/2/3 is the — )■ A isospin transition factor. Electromagnetic 
gauge invariance requires that (?/^r"^_j^^(pA, ?) = 0. The coupling constants Qi = QiiQ^ = 0) for 
i = 1,2,3 can be related to the magnetic, electric and Coulomb components of the 'jNA vertex by 
91 = dM, 92 = Qm + Qe, 5'3 = fi'c [3 [95]. The corresponding conjugate vertex can be obtained from the 
relation [3 [TOO] 

rA%7v(PA, q) = 70 [r"^_,^(pA, g)]So, (39) 

where pa and q are the momenta of the incoming A and outgoing photon. 

The amplitude of the box diagram with a A intermediate state can then be written as [TJ |95] 

-M(7^ = -ze^l^ Vi/rAF(gi,0)A^(g2,0), (40) 

where L^^ is the leptonic box tensor of Eq. (132|) . and the crossed-box is obtained by crossing symmetry, 
Eq. (125|) . The A hadronic tensor is given by 

= '«7v(p')rA%jv(P + ?i'-?2)^a/3(p + gi,MA)r^J^^^(p + gi,gi)M7v(p), (41) 
where the A propagator is given by 

Sa^{pA,M^) = -Sf{p^,M^)V%\pa), (42) 



r°^^^A(PA,g) 



2Ml V 3 



+ ^72(Q')[gX-^7"^g-PA] + 



Ma 



W (7X 
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where the projection operator 



Vjp (pa) = gap - -7a7/3 - ^ O^ATaPA/? + lfi^AV/\a) (43) 

ensures that only spin-3/2 components are present. Suppression of the unphysical spin- 1/2 contributions 
also leads to the condition on the vertex p/^oX'^N^Aiv^i l) — 0- Note that the mass difference Ma — M 
between the initial and intermediate hadronic states renders the integral (14T|) IR-finite, so that in 
contrast to Eq. (130|) . a photon mass parameter A is not needed to regulate divergences. 
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Figure 14: Sum of the nucleon and A(1232) contributions to the TPE correction to the ep 
elastic cross section for two values of the cut-off mass Aa- Figure taken from Ref. [7]. 

In Refs. [3 [95] the three 7A^A transition form factors gi{Q^) = Qi F^{Q^) (z = 1, 2, 3) were assumed 
to have dipole shapes, Fy{Q'^) = (1 + Q^/A^)"^, with a dipole mass Aa = 0.84 GeV for each. The 
electric and magnetic couplings were taken to have the values qm = 7 and qe = 9 [7], obtained from 
a K-matrix analysis of pion photoproduction data |99j. A more realistic vrA^ coupled channel quasi- 
potential study [101] gives similar values, qm = 6.3 and qe = 3.4. The Coulomb coupling gc is not as 
well constrained as the electric and magnetic couplings, and Kondratyuk et al. [7j considered the range 
between gc = —2 to 0. Tjon et al. [HS] used an estimate from the nucleon-A E2/M1 transition strength, 
gc = 5.8, but found the sensitivity to variations of gc small. Similar values for the couplings were also 
used by Zhou et al. |100] . although defined with different normalizations. The sensitivity to gc was 
found to be weak at low and e, becoming stronger at forward angles and large {Q^ > 3 GeV^, 

£-1) [inn]. 

The e dependence of the relative A and nucleon TPE contributions to the elastic cross section from 
Ref. |7], for gc = 0, is shown in Fig. [H] at = 1 and 3 GeV^. The most striking feature of A 
corrections is its positive sign, in agreement with the early estimates ^2} HU |l5l |?6J, which is opposite 
to that of the nucleon, and the corresponding negative slope in e. This has the effect of attenuating the 
(negative) nucleon TPE correction, making it somewhat smaller especially at backward angles. The 
magnitude of the A correction is considerably smaller than the nucleon, however, so that the general 
features of the TPE correction in Fig. [13] are not affected by the A. 

These features are also largely insensitive to the details of the 7A^A form factor, particularly at low 
Q^, as the comparison of the results in Fig. [TUwith cut-offs Aa = 0.84 and 0.68 GeV illustrates [7]. At 
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Figure 15: Reduced cross section, scaled by the dipole form factor squared, showing the 
effect of adding TPE corrections to the Born cross section. The intermediate state includes 
a nucleon and the indicated hadron resonances. The curves for = 2.64, 4 and 6 GeV^ 
have been shifted vertically by —0.04, +0.04 and +0.09, respectively, for clarity, and are 
compared with data from Refs. pTl [T3]. The nucleon-only result is slightly larger than the 
+ A(P33), lying close to the full calculation. Figure taken from Ref. |102] . 



higher Q^, on the other hand, both the magnitude of the A correction and its dependence on the form 
factor model increase, although for very large the pure resonance calculation itself becomes more 
questionable. 

4.4.2 Nucleon resonances 

If the A(1232) resonance makes a non-negligible contribution to the TPE correction, at least in some 
kinematic regions, the question naturally arises whether other, higher-mass resonances could also play 
some role. Kondratyuk and Blunden [102\ extended the formalism of Refs. O [7], generalizing it to 
include the full spectrum of the most important hadron resonances as intermediate states involving spin 
1/2 and 3/2 resonances. The masses of the resonances and their nucleon- photon coupling constants are 
based on dynamical multichannel calculations [991 1103[ 1104] of nucleon Compton scattering at low and 
intermediate energies. The resonance TPE effects turn out to be not too sensitive to the details of these 
models. 

In general the contributions of all the heavier resonances are much smaller than those of the nucleon 
and A (P33) [7]. However, there is an interesting interplay between the contributions of the spin 1/2 and 
spin 3/2 resonances, which is analogous to the partial cancellation of the two-photon exchange effects of 
the nucleon and A intermediate states, found in Ref. [7]. Notwithstanding the smallness of the resonance 
contributions, their inclusion in the TPE diagrams leads to a better agreement between the Rosenbluth 
and polarization transfer data analyses, especially at higher values of the momentum-transfer squared 

The total TPE correction is given by the sum of the separate hadron contributions, 

S = + 6^ + 6^'' + 6^'' + 5^" + 5^" + (44) 
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The coupling constants in the vertices are taken from the Dressed K-Matrix Model, whose essential 
ingredients are described in Ref. |99j. The calculated two-photon corrections to the reduced cross 
section are displayed in Fig. [T3 The one-photon exchange cross sections are shown by the dotted lines. 
The cross sections including additional TPE corrections are shown by the dashed lines for the sum of 
the nucleon and A contributions, and by the solid lines for the full result with all resonances. In general, 
each resonance two-photon correction is proportional to a sum of squares of the nucleon-photon coupling 
constants of that resonance. This sets the scale of the magnitude of the resonance contributions. 

As an example at the moderately high value = 4 GeV^, the TPE corrections from the included 
resonances can be classified by their signs and orders of magnitude as follows. For < e < 1, the 
corrections change smoothly between the values 

-5.0 < 5^ < 0%, 1.9 > 5^ > 0%, -0.7 < 5^^^ < 0%, 

-0.3 <, 6^'' <, 0%, -0.15 <, 5^" <, 0%, 0.06 ^ 5^" ^ 0%, 0.01 ^ 5^'' ^ 0%, 

listed in order of decreasing magnitude. Inclusion of the excited state resonance contributions therefore 
effectively reduces the nucleon elastic TPE correction by ~ 15% at this value. At lower the 
resonance contributions are even less important relative to the nucleon-only contribution. 

Figure [15] shows that at low to moderate the total TPE corrections as a function of e lie between 
those of the nucleon and A intermediate states. In addition to the dominant nucleon and A contribu- 
tions, the gives the most important correction among the remaining resonances. This is consistent 
with the well-known prominence of the Di^ in the second resonance region of the Compton scattering 
cross section {e.g. see Ref. |103[ 1104] and references therein). 

Inclusion of contributions of intermediate states with masses larger than ~ 2 GeV becomes imprac- 
tical within a hadronic approach when one moves into the deep-inelastic continuum. Here it becomes 
more efficient to use either partonic degrees of freedom, discussed in Sec. 14.61 or dispersion relations, 
discussed in Sec. 14.71 

4.5 Scattering from a point target 

It is useful to examine the TPE corrections to elastic scattering from a structureless pointlike target, 
such as a /i"*". or in the limit of hard scattering from quarks in the nucleon. The pointlike limit is 
realized by setting Fi = 1 and F2 = in the electromagnetic current operator appearing in Eq. fl33]) . or 
equivalently by replacing Tf^{q) — t- 7^^. The reduced cross section is then = t + e. In this case, we 
take the IR divergent part of the TPE correction to be the model-independent Maximo n-Tj on result of 
Eq. fl2B]) : the remainder is denoted as ^hard, viz. 5^^ = 5iR(MTj) + ^hard- 
Figure [12] shows the result for ^hard as a function of e for various values of ranging from to 00. 
Clearly 5hard is a positive quantity which generally increases with Q'^ between the two limiting cases. 
In both limits simple expressions for 5hard exist in terms of a single variable. For a pointlike target the 
results should exhibit a scale invariance under the appropriate dimensionless variable, which can be 
taken to be r and e. 

(i) The limit r — t- is the same as either — t- or the static limit M — )■ 00 with finite. In 
this case 5hard is a function of a single variable. Making use of the asymptotic expansions of the 
Passarino-Veltman functions, a particularly simple expression exists in this limit, namely 

^hard = —, (45) 

X -I- 1 

with X = {1 + e)/{l — e). This static limit is also realized in the second Born approximation, 
which we discuss below. It agrees with the expression given in Ref. [39], provided one identifies 
sin(6'/2) = 1/x in this limit. 
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(ii) The limit r — )■ cxd is the same as either — )■ cxo or the hmit M — > with finite. Once again 
(5hard is a function of a single variable only, and can be expressed as 
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hard 



vr x2 + 1 
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x + 1 
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+ X 



vr^ + In^ 



X + 1 



+ ln2 



X 



In 



(46) 



This expression, plus the soft IR contribution 5iR(MTj), agrees with that for scattering from 
massless quarks given in Eq. (27) of Ref. [86], and also in earlier work. However, those results 
were not expressed in terms of a single variable. 



The behavior of TPE for hard scattering from a pointlike target stands in stark contrast to the 
ep elastic scattering results, shown in Fig. [131 which is larger and has the opposite sign at high . 
However, as one decreases Q^, the ep TPE correction becomes smaller, eventually switching sign from 
negative to positive as decreases. For — )■ 0, the hadronic result reproduces the pointlike limit 
shown in Fig. [16] (and given in Eq. (l45l) ) - independent of the anomalous magnetic form factor F2. (Note 
that (5iR(MTj) — !■ 5ir(MoT) in the limit — 0, so that 5 — )■ (5hard)- This suggests that the magnetic 
interaction is unimportant in this limit. This is consistent with the low-Q^ analysis of Ref. [6^ . 

Kuraev and Tomasi-Gustafsson [105j have suggested that e/i scattering should be a function of x, and 
that it represents an "upper limit" for the hadronic TPE corrections. In the context of the discussion 
here, while this may be a positive upper bound, at high the full calculation gives a correction larger 
in magnitude but opposite in sign. 

To fully understand the underlying reasons for this behavior, one can invoke the second Born ap- 
proximation. In the static limit, equivalent to the case where the proton has infinite mass, the scattering 
amplitude for TPE coincides with the amplitude for electron scattering in a Coulomb potential in second 
Born approximation. This limit was first considered by Dalitz |106j . 

In Fig. [17] the TPE correction is shown as a function of Q for a typical electron energy involving a 
range of low values |107j . The TPE contribution is compared to the contribution involving only the 
Coulomb distortion of the electron, arising from a second soft photon, calculated according to Ref. |108] 
in second Born approximation. The same contribution, but for a point nucleus, is given for comparison. 
Remarkably good agreement between the two approaches is obtained at forward angles. Figure [17] 
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Figure 17: Relative contribution of TPE to elastic ep scattering at E = 160 MeV. The results 
in second Born approximation account for the Coulomb distortion (exchange of second soft 
photon) only. Figure taken from Ref. [107j . 



shows that at forward angles the TPE contribution is dominated by the Coulomb distortion, while at 
backward angles the exchange of two hard photons contributes appreciably. This is consistent with the 
observations about the pointlike limit, and the independence of this limit from F2. 

The behavior shown in Fig. [13] (left) can also be understood from the second Born approximation 
[109] . The second interaction provides a "focusing effect" — accelerating electrons towards the target. 
For a pointlike target this should increase scattering at backward angles. However, because of an 
increased momentum transfer there is a competing effect from a reduction in cross section due to the 
proton electric and magnetic form factors. At larger this reduction wins out, and the total cross 
section is reduced at backward angles. For positrons, the opposite effect is expected. 

4.6 High partonic models 

In the regime of high Q^, two approaches to the TPE effect on elastic scattering have been taken by 
different groups. In Refs. [551 IBS], the hard scattering part of TPE was studied in a partonic approach 
using different models for generalized parton distributions (GPDs). In this approach it is assumed that 
both photons interact with the same quark. The results of these calculations on the cross sections and 
other observables have been examined in detail in Ref. |110j . and the details will not be duplicated here. 



(a) 1 (b) (c) i 

Figure 18: Typical pQCD diagrams for (a) one-photon exchange, (b) leading order TPE 

involving 1 hard gluon, and (c) subleading order TPE involving 2 hard gluons. Figure 
adapted from Ref. |lllj . 



Our observations about the pointlike limit in Sec. 14.51 suggest that hard TPE corrections involving 
both photons interacting with the same particle have the opposite sign to the hadronic calculations 
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at backward angles. Borisyuk and Kobushkin [lllj studied TPE in the framework of perturbative 
QCD (pQCD). In this approach, it turns out that the most important diagrams for backward angle 
scattering are the ones where the two photons interact with different quarks. This allows the possibility 
to generate a hard TPE correction that is negative at backward angles. The argument given by Borisyuk 
and Kobushkin is as follows. 

Referring to the diagrams in Fig. [181 one-photon exchange diagrams like (a) need 2 hard gluons to 
turn the momentum of all 3 quarks, and are therefore of order aa1/Q^ , where as is the strong coupling. 
The TPE diagrams like (b), involving 2 quarks, need only 1 hard gluon, and are of order a^a^jQ^ . 
Therefore the ratio of two-photon to one-photon exchange is not of order a, as one might expect, but of 
order a/us- By contrast, diagrams like (c), where both photons interact with the same quark, require 
2 hard gluons, and are therefore suppressed by an additional factor of a^. 
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Figure 19: TPE amplitude vs. for two values of e. The dashed curves are the hadronic 
calculation using two different form factors, and solid curves are the pQCD results. Figure 
taken from Ref. |111] . 



Numerically, Borisyuk and Kobushkin |111] express the cross section correction in terms of the 
relative TPE amplitude 6Gm/Gm, where Gm = G^°^'^ + SGm, which has linear e dependence. This 
amplitude grows logarithmically with Q^, reaching 3.5% at = 30 GeV^. This may offer an interesting 
avenue to merge the hadronic approach, valid at low to moderate Q^, with the pQCD approach at 
high Q^. Figure [19] shows how the hadronic and pQCD results may connect approximately at around 

= 3 GeV^ for two values of e. A complementary pQCD analysis by Kivel and Vanderhaeghen |112] 
also finds a comparable effect. 

Two-photon effects also manifest in the production of lepton-antilepton pairs from a hadronic target 
by real photons. The interference of one- and two-photon exchange amplitudes leads to a charge 
asymmetry term that can be calculated explicitly in the large-t limit. This was studied using pQCD 
by Hoodbhoy |113] in the region of high center of mass energy s (s ^ — t ^ M^). Typical diagrams 
involving 2 photons and 1 hard gluon exchange are shown in Fig. analogous to the corresponding 
diagrams in elastic ep scattering discussed previously. The lepton pair asymmetry for select angles 
in the lepton pair center of mass is also shown as a function of p = —t/M'^. Note that Hoodbhoy's 
calculation and observations about the leading contribution involving different quarks actually predates 
the work of Refs. [TTT|[TT2]. 
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Figure 20: Typical TPE diagrams for lepton pair production involving two photons and one 
hard gluon (left). The graph shows the resulting lepton pair asymmetry at select angles 
(right). Figures taken from Ref. |113j . 



4.7 Dispersion relations 

One of the sources of model dependence in the calculation of the real parts of the TBE amplitudes arises 
from the intermediate state in the box diagrams being off-shell. For the nucleon elastic contribution, for 
example, this can be parametrized through the off-shell dependence of the nucleon electromagnetic form 
factors, where one of the nucleons is off its mass shell, = {p + qi^ M^, see Fig.HOl In this case the 
vertex form factors are in general functions of two variables, qf and W'^, which in the limit W"^ — )■ 
reduce to the on-shell values. Outside of a field-theoretical calculation of the amplitudes from first 
principles, however, it is not possible to determine the (unphysical) dependence on independently 
of a specific model. This model dependence invariably introduces additional uncertainties into the 
calculation. 

One way to avoid the off-shell ambiguities is through the use of dispersion relations |114j . which 
are based on general principles such as unitarity and analyticity of scattering amplitudes, and can in 
principle relate real and imaginary parts of amplitudes for on-shell processes. The conservation of total 
probability in a scattering process implies that the scattering matrix, S, is unitary, S^S = 1. In general 
the S-matrix element Sfi = {f\S\i) between initial state i (with total momentum pi) and final state / 
(total momentum pf) can be written in terms of the invariant amplitude as 

Sf, = 6f, + t{27cf6\pf-p,)Mf,. (47) 

Unitarity then requires that 

2ZmMfi = J2 I dpM*^fMni, (48) 

where dp is the covariant phase-space factor for a complete set of intermediate states n. In the limit 
of forward elastic scattering, i = f, the product of amplitudes in (HSj) becomes diagonal, and can be 
expressed as a directly measurable cross section. This allows loop integrations to be evaluated in terms 
of on-shell states, so that empirical data can be used as input in the calculation of the imaginary part. 

The Cauchy integral formula allows the amplitude Aifi{s) to be determined as a function of a real 
s-channel variable s from the structure of its singularities in the complex plane. The real and imaginary 
parts of Aifi{s) can be related by 

neM,M ^ r ,s'^-^^^P^, (49) 

TT J-oo S' - S 

where V denotes the principal value of the integral. 
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The dispersion approach has been apphed in various studies in the hterature of TPE corrections to 
elastic ep scattering, starting with some early attempts in the late 1950s [3111121113] to compute 0{e^) 
corrections to elastic ep cross sections. More recently, Borisyuk and Kobushkin |115] used dispersion re- 
lations to compute the contribution to TPE amplitudes arising from nucleon elastic intermediate states, 
in terms of elastic on-shell nucleon form factors parametrized as a sum of monopoles [5]. Interestingly, 
the results of the dispersive calculations were found to be very similar numerically to the direct loop 
computation [6l 1116] . especially at low Q^. This suggests that the prescription of taking on-shell values 
for the half off-shell form factors in the direct calculations provided a very good approximation. At 
larger (Q^ > 6 GeV^) the agreement between the two approaches deteriorates, possibly hinting at 
the need to account for the off-shell dependence of the form factors; at such values the reliability 
of the hadronic approach itself may be more questionable. In Ref. |117] the authors applied the same 
method to also compute the TPE corrections to the pion form factor (see Sec. 16.41 below), again finding 
very similar results to the direct loop calculations 

Elastic scattering from hadrons at forward angles was analyzed some time ago within a dispersive 
approach by Bernabeu and collaborators |118[ I119[ I120[ 1121] , and more recently by Gorchtein |122] . We 
sketch here the relevant details (in the notation of this article), noting the overlap with the discussion 
of forward-angle parity-violating 7Z interference in Sec. 17.31 In the forward angle limit {qi ~ —q2 in 
Fig. [TU]) the imaginary part of the TPE amplitude, by virtue of the intermediate state being on-shell, is 
related via the optical theorem fHS]) to the inclusive electromagnetic structure functions of the nucleon. 
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where I = k — qi is the momentum of intermediate state electron, Ei = ^J{k — qiY + ml its energy, 
and Q\ = —qf is the virtuality of exchanged photon. The factor AirM arises from the definition of 
the hadronic tensor. The leptonic tensor in (150|) is given by L^^, = Ueik')'-y^{J + 771^)^^ Ue{k) , where 
the scattered lepton momentum has been set to k' = k — q, with q kept small but finite. The nucleon 
hadronic tensor can be written in terms of the inclusive Fi and F2 structure functions as 
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dropping terms proportional to q^^v which vanish when contracted with L^y. Using Eq. ( l35|l the relative 



correction from TPE to the elastic cross section can then be written as 
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where the upper limit on the Ql integration is Q^^^ax = 2M_E'(1 — W"^ /s), with s = + 2ME the 
total center of mass energy squared, and t = q^. Note that the factor t in Eq. (!52|) implies that formally 
this contribution vanishes at t = 0. 

The real part of 5^^{E) at a given electron energy E is then evaluated from the imaginary part 
using the fixed-t dispersion relation 



ne5yy{E) 



TT 



dE 



, Xm 5-f^{E') 



E' -E 



(53) 



Note that the contribution to the integral from the region E' < corresponds to the crossed-box 
diagram in Fig. [TOl Invoking the crossing symmetry property discussed earlier in Eq. f l25|) . one can 
write 

2E f°° 1 

neS,,{E) = —V dE'^j^-^Xm6,,{E'). (54) 
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Bernabeu et al. [llSj I119[ I120[ 1121] applied this formalism to the case of elastic electron and muon 
scattering from the spin-0 ^He and ^^C nuclei at low energy, focusing in particular on the continuum 
contributions. The electron-positron charge asymmetry for ^He, which depends on the interference of 
the one- and two-photon exchange amplitudes, was found to be of the order of 1%, and about 4-5 smaller 
for ^'^ |120j . The effects of the finite size of the nucleus were found to be small for kinetic energies up 
to ~ 50 MeV, but increased to ~ 30% of the total forward amplitude at energies ~ 100 MeV |121] , and 
somewhat more significant for muon scattering. The results of these studies were subsequently used to 
constrain the subtraction constant in the dispersion relation for the two-photon exchange amplitude at 
nonzero values of the momentum transfer squared t. 

More recently, Gorchtein |122] applied the forward dispersion formalism to study the effect of TPE 
on the e^p to e~p cross section ratio at higher energies. Using phenomenological structure functions 
in the deep-inelastic scattering region |123] . the ratio of e+p to e~p cross sections is found to be larger 
than unity, with effects ranging from <, 1.5% for incident energy E = ?> GeV to ^ 2.5% for £^ f» 10 GeV. 

Overall, the dispersion approach provides a unique method to compute amplitudes from empirical 
inputs, in principle free of model dependent ambiguities associated with off-shell effects of intermediate 
states, particular at high mass W"^. On the other hand, this virtue is at times somewhat negated 
by the need to know the inputs at all energies, which are often not well determined at very high 
energy. Uncertainties in the high energy behavior of the input cross sections can in some cases be as 
large as the off-shell uncertainties in the direct, non- dispersive approaches. Furthermore, at scattering 
angles away from the forward limit the imaginary parts of TPE amplitudes are no longer related to 
inclusive structure functions, but are given by nonforward virtual Compton scattering amplitudes, for 
which there are considerably fewer phenomenological constraints. The dispersive framework therefore 
provides a useful complement to the other methods discussed in this section. 

5 Impact of two-photon exchange on observables 

Having outlined the theoretical developments in the computation of TPE corrections to elastic scatter- 
ing, in this section we review the implications of the corrections for observables. Starting with the most 
topical case of TPE effects on the proton electric to magnetic form factor ratio (Sec. 15.11) . we also con- 
sider the impact that TPE corrections have on electron-proton to positron-proton elastic cross sections 
(Sec. 15.21) . on the global analysis of form factor data (Sec. 15.41) . as well as on polarization observables 
(Sec. 15.31 and 15.51) . To simplify the discussion, unless otherwise indicated, we restrict the hadronic 
calculations to the nucleon elastic intermediate states only (which provides a reasonable approximation 
to the total hadronic contribution, see Sec. I4.4.2p . 

5.1 Cross sections and Rosenbluth separations 

The effect of the TPE corrections on the Rosenbluth separations, calculated within the hadronic for- 
malism of Sec. 14.31 is illustrated in Fig. |2T1 The data points are the reduced cross sections an, scaled 
by the square of the dipole form factor (^71) multiplied by r, from SLAG [M] and the Jefferson Lab 
Super- Rosenbluth experiment [H]. Here the Born- level prediction (dotted line) is obtained using the 
form factor parametrization of Ref. [97] in which Ge is fitted to the polarization transfer data, while 
the full result (solid) includes the TPE corrections. The TPE contributions yield a significant increase 
of the slope, with some nonlinearity evident at small e. The corrected results are clearly in better 
agreement with the data, although some residual difference between the Rosenbluth and polarization 
measurements remain at the higher values. 

To obtain a rough estimate of the influence of these corrections on the electric to magnetic proton 
form factor ratio R = ^pGe/Gm, Blunden et al. [5l |6] considered a simplified analysis in which the 
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Figure 21: Comparison of the Born (dotted) and TPE (solid) calculations to measurements 
of aR/{TGl) for several values of Q"^: (a) SLAG data [34j at = 3 25 - 6 GeV^; (b) 
Jefferson Lab data [14J at = 2.64 — 4.1 GeV^. The Born cross section is evaluated using 
the form factors from Ref. [HZ] which use PT extractions of G^;. The curves in (b) have been 
shifted by (+1.0%, +2.1%, +3.0%) for = (2.64,3.2,4.1) GeV^. Figure adapted from 
Ref. ig. 



effective e slope was approximated by a linear fit to the full TPE corrections. Such an analysis may 
not provide an accurate estimate at very low e or high where strong deviations from linearity arise; 
however, it is still instructive to obtain an estimate of the effect on R by taking the slope over several 
ranges of e. Blunden et al. [6] fit the relative correction to the Born cross section (1 + 5) to a linear 
function of e, of the form a + 6e, for each value of at which the ratio R is measured. The corrected 
reduced cross section in Eq. (ITU]) then becomes 



where 



1 + ^ (i?2 [1 + eh/ a] + 



(55) 



1 + eb/a ^ ' 

is the "true" form factor ratio, corrected for TPE effects, and R is the "effective" ratio, contaminated by 
TPE. Note that in Eqs. flS^ and fISB]) the term quadratic in e has effectively been linearized by averaging 
e over the range fitted, e ^ e. The approximation in Eq. fISS]) is reasonable provided 6/a ^ 1, which is 
more valid at high values. 

Gonsidering two ranges for e, namely a large range e = 0.2 — 0.9 and a more restricted range 
e = 0.5 — 0.8, the resulting shift in R is shown in Fig. [22l together with the polarization transfer 
data. The linear e approximation to the TPE correction should be better for the latter, even though 
in practice experiments typically sample values of e near its extrema. 

The effect of TPE on R is clearly significant. In particular, the corrections have the proper sign 
and magnitude to resolve a large part of the discrepancy between the LT and PT techniques. While 
the early analysis of Ref. using simple monopole form factors found a shift similar to that in for the 
e = 0.5 — 0.8 range in Fig. which resolves around 1/2 of the discrepancy, the nonlinearity at small e 
makes the effective slope somewhat larger if the e range is taken between 0.2 and 0.9. The magnitude of 
the effect in this case is sufficient to bring the LT and PT points almost to agreement. A more detailed 
analysis, where the TPE correction is applied directly to the experimentally measured cross sections, is 
discussed in Sec. I5.4.2I 
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Figure 22: The ratio of proton electric and magnetic form factors iipGe/Gm measured using 
LT separation (open diamonds) [12] and polarization transfer (PT) (open circles) [TB]. The 
LT points corrected for TPE are shown assuming a linear slope for e = 0.2 — 0.9 (filled 
squares) and e = 0.5 — 0.8 (filled circles) (offset for clarity). Figure adapted from Ref. [B]. 



Note that for the high points, the value of iipGe/Gm decreases but the uncertainty on the ratio 
increases when TPE corrections are applied. This is because the quantity directly constrained by the 
experiment is the slope of the reduced cross section, which is related to {(IpGe/GmY, rather than 
(J^pGe/Gm- So while the absolute uncertainty on the slope, and thus {^ipG e / G m)'^ , is unchanged, the 
uncertainty that this translates to for UpGe/Gm depends on the value, and increases when UpGe/Gm 
decreases. 

The results in Fig. |22] are representative of most calculations of the TPE corrections within the 
conventional hadronic approach. Qualitatively similar results are also obtained in the GPD-based 
approach [53 ES], whose applicability is limited, however, at low and e values (c/. Fig. (TH]). Note 
that since the hard electron-quark scattering amplitude gives a negative slope in £, Fig. Uni agreement 
with data suggests an important role played by soft physics. 

5.2 e^p/e~p ratios 

As previewed in Sec. 13.3. direct experimental evidence for TPE effects can be obtained by comparing 
elastic e"^p and e~p cross sections. The ratio of these is defined as 

^e+p |A^f |2 + 27^ef^^f*^^!;;) 

a^-p |^1^-|2 + 27^e(^^^'*^<^-) ' ^ ' 

where a^^^ = da'^^'^ / dVt. Whereas the electron Born amplitude M.^^ changes sign under the interchange 
6^—7-6+, the TPE amplitude ■M'^^ does not. The interference of the Born and TPE amplitudes therefore 
has the opposite sign for electron and positron scattering, so that the cross section ratio can be written 

^ 1-25, (58) 

where 5 is defined in Eq. ( !36|) . The hadronic TPE calculation of 5 is illustrated for a range of values 
in Fig. [131 Since the finite part of the TPE contribution is negative over most of the range of e, one 
would expect to see an enhancement of the ratio of e"*" to e~ cross sections. 
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Existing data on elastic e~p and e~^p scattering are sparse, although some constraints exist from early 
measurements at SLAC [SSI El], Cornell [TOj, DESY [Tlj and Orsay [73j, as discussed in Sec. I3.3.I1 The 
data are predominantly at low and at forward scattering angles, corresponding to large e {e 0.7), 
where the TPE contributions to the cross section are small 1%). Nevertheless, the overall trend in 
the data reveals a small enhancement in W^^^ at the lower e values, as shown in Fig. [71 

The measurements of R'^^'^ are compared to the hadronic TPE calculations [6] in Fig. [231 which 
generally predict a slight enhancement in the ratio at the experimental kinematics. Overall, the TPE 
corrections agree reasonably well with the data, although the experimental uncertainties are quite large, 
especially where there are indications of a nonzero TPE effect. Interestingly, the GPD-based calculations 
[551 IBS] predict a suppression of i?^^*^ at large e in the ~ few GeV^ range. Better quality data, 
particularly at backward angles, where an enhancement of up to ~ 10% is predicted, are needed for a 
more definitive test of the TPE mechanism. 
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Figure 23: Measured ratio of elastic e~^p to e~p cross sections (crosses), compared to the 
hadronic TPE calculations [6] (diamonds), presented as a function of e (left) and (right). 



In the last few years, three new experiments have been initiated to compare e~^p and e~p scattering. 
The first took data in 2009 at the VEPP-3 ring in Novosibirsk [77], which uses an internal target in a 
positron/electron storage ring to extract the ratio at = 1.6 GeV^ and e ~ 0.4 [124] . The analysis 
yields a raw e~^p/e~p ratio of 1.056 ± 0.011, but corrections for charge-dependent bremsstrahlung, 
estimated to be ~ 3%, have to be applied [125j . The completed experiment is statistics limited, with 
estimated systematic uncertainties of 0.3% [125] . and plans for a follow-up measurement with higher 
statistics at similar or somewhat lower are being evaluated. 

Another experiment [78], using a mixed beam of e"*" and e~ generated via pair production from a 
secondary photon beam, recently completed data taking at Jefferson Lab. In this experiment, e~p and 
e~^p elastic scattering are measured simultaneously, using detection of both the scattered lepton and 
struck proton to reconstruct the initial lepton energy. Cross sections can be measured for 0.5 < < 
2.0 GeV^ and e ^ 0.2, and the analysis is currently underway. The use of a lepton beam with a large 
range of energies allows for a mapping out of the e dependence at fixed values over the kinematic 
range of the experiment. 

Finally, the OLYMPUS experiment p9] will use the DORIS lepton storage ring at DESY to make 
measurements with a fixed lepton energy, yielding e^p/e~p ratios at several points from = 2.2 GeV^, 
e = 0.35, to = 0.6 GeV^, e = 0.9. The experiment is currently being installed and is scheduled to 
run in 2012. 
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5.3 Polarization measurements 



While the results of the modern TPE calculations give a clear indication of a sizable correction to the 
LT-separated data, the obvious question which arises is to what extent does TPE affect the polarization 
transfer results themselves. The expectation is that since the FT measurements involve ratios of cross 
sections, most of the radiative effects, including TPE, should cancel. 

The polarization transfer experiment involves the scattering of longitudinally polarized electrons 
from an unpolarized proton target, with the detection of the polarization of the recoil proton, e + p — )■ 
e + p. (The analogous process whereby a polarized electron scatters elastically from a polarized proton 
leaving an unpolarized final state gives rise to essentially the same information.) The calculation of 
the TPE corrections to the cross sections with longitudinally or transversely polarized recoil protons 
follows that of Sec. 14.31 with the spin traces evaluated using the explicit expression for the spin- vectors 
of the incident electron and recoil proton in Eqs. ( |T2l) and ( |T3i) . 

In the standard radiative corrections using the Mo-Tsai prescription [371 EH], the corrections to 
the transverse and longitudinal polarization are identical, so that no additional corrections beyond 
hard bremsstrahlung are necessary [15]. For the TPE corrections, in analogy with the unpolarized case, 
Eq. fl35|) . the spin- dependent TPE corrections to the longitudinal (Sl) and transverse {6t) cross sections 
are defined as the finite parts of the TPE contributions relative to the IR expression from Mo and Tsai 
[2Zl EH] in Eq- ( 12^ . which are independent of polarization, 

5l,t = h,T - <5ir(MoT). (59) 

In the polarization transfer measurements of the form factor ratio (IpGe/Gm, one usually measures 
the ratio of the longitudinal or transverse cross section to the unpolarized cross section, Pl or P^, 
respectively. The TPE correction to the polarization transfer ratio can therefore be incorporated as [6] 

Pl,t 1 + 5l,t 

= 1 , T ' (60) 

where -P^t polarized Born cross sections, and the correction to the unpolarized cross section 6 is 

given in Eq. ( 136|) . Polarized target measurements make similar comparisons of beam-target cross section 
asymmetries for different target spin orientations, yielding two different combinations of longitudinal 
and transverse spin asymmetries. 

Taking R = —^p^Jt{1 + e)/2e (Pt/Pl) to be the corrected ("true") electric to magnetic form factor 
ratio (see Eq. f l56|) ). the experimentally measured PT ratio is 

Inverting Eq. ( 16T|) , the shift in R is illustrated in Fig. [2U relative to the uncorrected results and the LT 
separated data. Clearly the effect of TPE on the form factor ratio is a very small, ^ 3% suppression at 
the larger values, which is well within the experimental uncertainties. Note that the shift in R in 
Eq. I ET\i does not include corrections due to hard photon bremsstrahlung, which are part of the standard 
radiative corrections. Since these would increase both the numerator and denominator in Eq. (1611) . the 
correction in Fig. |24] represents an upper limit on the shift in R. 

Although the TPE effects on the PT ratio R are very small, this is mostly because those mea- 
surements were typically performed at large e {e ^ 0.7 — 0.8) (TB], where the TPE corrections are 
minimal. At lower e, or for e ~ 1, the effects on the polarizations can be significant. This can be seen 
in Fig. [25], which shows the longitudinal and transverse polarizations from nucleon elastic intermediate 
states relative to the Born terms. The correction Sl to the longitudinal cross section is approximately 
the same as the correction 5 to the unpolarized cross section, so that TPE correction to the longitudinal 
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Figure 24: Proton electric to magnetic form factor ratio obtained from the PT measurements 
[16j, with (sohd circles) and without (open circles) TPE corrections. The corrected values 
have been offset for clarity. The LT-separated ratio (open diamonds) from Fig. |22] is shown 
for comparison. Figure adapted from Ref. 



polarization is extremely small. In fact, at extreme backward angles {e = 0) the unpolarized and 
longitudinal corrections 6 and 6l must be identical [6j, and there is no TPE contribution to P^. In 
contrast, the combined effects of an increasing St at large and S becoming more negative mean that 
the correction to the transverse polarization P^ is enhanced at backward angles, and grows with Q^, as 
Fig. [2S1 illustrates. 

In a recent Jefferson Lab experiment the ratios R and Pl/P^ were measured at three e values for 
a fixed = 2.5 GeV^ [67j. In the Born approximation both of these should be independent of e, 
so that any observed angular dependence would be an indication of TPE contributions. The results 
for the longitudinal polarization ratio in Fig. [25] (left) are consistent with no TPE effect at e = 0.635, 
but indicate an ^ 2% enhancement at the larger e = 0.785 value (note that an overall normalization 
uncertainty of 0.01 applies to the data). While a rise as £ —> 1 is predicted at larger Q^, the magnitude 
of the effect is difficult to reconcile with the calculation [6j and possibly suggests the need to include 
higher-mass intermediate states in the TPE amplitude. The results in the GPD-based approach [55| [86] 
give a similarly small effect, but with opposite sign compared with the hadronic calculation |6]. 

Data taken on the PT ratio R ai = 2.5 GeV^ are also consistent with no significant e dependence 
over the range 0.15 ^ e ^ 0.8, and hence in the transverse ratio Pt/Pt- The hadronic calculation with 
elastic nucleon intermediate states predicts an ?a 2 — 3% enhancement, which once again may indicate 
a role played by the higher-mass resonances; the A(1232) contribution, for example, may cancel some 
of the rise at e — )■ (Sec. I4.4.ip . It is interesting to observe that the GPD-based model predicts a large 
decrease in the ratio R (or transverse polarization) at low e, opposite to the hadronic model. Further 
measurements of the e dependence of the Pt/P^ ratio over a range of values would be very helpful 
in constraining the magnitude and sign of the TPE effects. 

Finally, two recent analyses [65] [66] have attempted to extract from the new data individual contri- 
butions to TPE amplitudes. Because data exist for a limited set of values, and because the effect 
on the cross section is only inferred from the discrepancy between LT and PT measurements, such 
analyses at present require assumptions about the e dependence of the amplitudes. Although different 
decompositions of the amplitudes are used in Refs. [651 ES], the extracted amplitudes are found to be 
similar when the different definitions are accounted for. 
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Figure 25: Ratio of the finite part (with respect to the Mo-Tsai IR contribution ( 129|) ) of 
the Born+TPE correction relative to the Born term, for (left) longitudinal and (right) 
transverse recoil proton polarization, at = 1 (dotted), 3 (dashed) and 6 GeV^ (solid). 
The longitudinal polarization data are from Ref. [67], and have an overall normalization 
uncertainty of 0.01. Figure adapted from Ref. [6]. 



5.4 Global form factor analysis 

With the realization that the Rosenbluth separations of elastic cross section data were significantly 
affected by TPE, came the need to ensure that for reliable extractions of the nucleon form factors the 
TPE corrections must be incorporated in the analysis. Initial efforts focused on minimizing the impact of 
TPE in the extractions, based on the assumption that TPE was the likely explanation for the LT/PT 
form factor discrepancy, and that all observables (cross sections and polarization components) were 
affected at a similar level. Given that there are now calculations of TPE corrections that can explain 
the discrepancy and are consistent with all experimental constraints on TPE contributions, more recent 
efforts have focused on correcting the data using calculated TPE contributions in the combined analysis 
of cross section and polarization measurements. 



5.4.1 Phenomenological analyses of Rosenbluth and polarization data 

Because the Rosenbluth extractions at high are very sensitive to any small, angle dependent system- 
atic effects, it was generally believed that the polarization measurements provided a more reliable way 
to extract G^; at large Q^, where it has minimal contribution to the cross section. Section l3^ discussed 
various attempts to use the discrepancy to extract the TPE contributions to the data. However, when 
combining these measurements to extract the form factors, rather than the TPE, the exact result is 
sensitive to the way these inconsistent measurements are combined. 

Many early attempts to extract form factors using both unpolarized cross section and polarization 
transfer measurements simply assumed that Gm could be reliably extracted from the Rosenbluth re- 
sults, and Ge could best be extracted using polarization measurements of Ge/Gm, combined with the 
Rosenbluth extractions of Gm- This corresponds to the assumption, illustrated in Fig. [5l that whatever 
was deficient in the cross section measurements led to a linear e dependence which had no impact at 
e = 0, and thus the value of Gm extracted from Rosenbluth separations was unmodified. However, if 
the discrepancy is caused by TPE contributions, then it is e = 1 where these contributions must be 
negligible, and thus there will be both a large change in Ge and a smaller change in Gm- While the 
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size of the discrepancy suggests that the impact on Gm would be on the order of 2 — 3%, this is large 
compared to the precision with which Gm extractions were quoted and thus needs to be taken into 
account. 

In practice, different combined analyses made assumptions that led to somewhat different results. 
Because the polarization measurements were only available with precision at high Q^, some analy- 
ses [126] used Gm from LT separations, while taking Ge from LT separations at low and polarization 
measurements at high Q^. A slightly different approach was taken in Ref. [97j. Rather than fitting the 
extracted values oi Ge and Gm from a mix of cross section and polarization measurements, the slope 
of the reduced cross section was fixed from polarization measurements, and Gm varied to best match 
the reduced cross section over the entire e range. This corresponds to taking the slope predicted by 
polarization transfer in Fig. [5] and varying the normalization so that it gives the best reproduction, on 
average, of the cross section measurement. In this case, the slope will match the data at a weighted 
average value of e, typically reproducing the cross section at a point in 0.5 < e < 0.8. This is closer to 
what one expects from TPE, where the agreement would occur at e = 1, but still yields a value of Gm 
that is too small. 

One issue of concern is the fact that the TPE-corrected form factors will no longer provide a good 
parametrization of the form factors unless TPE contributions are explicitly included. This scenario can 
be seen from Fig. El where Gm is determined from the extrapolation of the cross section measurements 
to e = 0, and then the polarization measurements used to give the e dependence of the reduced cross 
section. The form factors thus obtained, if treated in the Born approximation, will yield the correct 
cross section only for scattering at 180°, and produce the maximum deviation from the measured cross 
section at small angles, where a large fraction of the high-precision measurements is taken. If these 
form factors are used as input in the analysis of other measurements, such as when calculating the ep 
elastic cross section in the analysis of quasi-elastic proton knockout from a nucleus [36l [STj 1127] and the 
deuteron elastic form factors |128j . then this error in the cross section can yield incorrect results. 

Note that this is not related to assumptions or errors in the extraction. The fact that TPE contri- 
butions are larger than previously thought means that it is not sufficient to take the Born form factors 
when one needs precise knowledge of the elastic cross section; the additional TPE contributions must be 
included along with the traditional radiative corrections to obtain the observed cross section. Because 
different analyses use different approaches to correct for the TPE contributions, care must be taken to 
apply consistent TPE corrections when extracting the form factors and obtaining the cross sections from 
these. One way to avoid these issues is to provide both an extraction of the TPE-corrected form factors 
and a fit to the total (Born+TPE) elastic cross section [36l 11291 1130] . This provides a parametrization 
of the cross section that includes the observed TPE contributions. 

5.4.2 Extraction of form factors including TPE 

Given that modern TPE calculations appear able to resolve most or all of the form factor discrepancy, 
the next step is to apply these corrections to the experimental observables and see if the results are 
consistent. If so, a combined extraction of the form factors from all available data can be performed. 

Such an analysis was performed in Ref. [I29J using the hadronic TPE calculations with nucleon elastic 
intermediate states [6], as outlined in Sec. 14. 3[ The hadronic TPE corrections have the advantage that 
they are expected to be reliable for all e values at low Q^, which is the kinematic region where many of 
the extremely high precision measurements have been taken. At higher values, the calculations is 
expected to be less complete, and it is important to test the corrections against the observed discrepancy. 
For ^ 2 — 3 GeV^, the TPE corrections to the cross section bring the form factor ratio extracted from 
Rosenbluth separations into excellent agreement with the polarization transfer measurements (middle 
panel of Fig. [26|l . At higher there is a small residual, systematic disagreement. Here a small 
additional phenomenological correction was applied, linear in e and with magnitude comparable with 
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estimates [71 [55], 1102] of higher-mass intermediate state contributions (right panel of Fig. [26 
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Figure 26: Comparison of polarization measurements (filled diamonds) and LT separations 
(open circles) with no TPE corrections (left), TPE corrections from Ref. |6j (center), and 
with the additional high-Q^ correction applied in Ref. [129\ (right). 



For the combined analysis of cross section and polarization measurements, TPE corrections were 
applied to the extracted cross sections. Most experiments assumed a 1-1.5% uncertainty due to radiative 
corrections, with the dominant contribution coming from TPE. Clearly, this was an underestimate of the 
uncertainty when no TPE corrections were applied, and was taken to be an appropriate uncertainty after 
applying the hadronic correction of the TPE effects. For the additional TPE contribution associated 
with higher-mass intermediate states, 100% of the correction was applied as an additional uncertainty to 
the cross section to reflect the impact of the poorly constrained TPE corrections at high Q^. While the 
TPE calculation [6J provides predictions for the impact on the polarization transfer measurement, most 
of the data are at large e, where the impact is extremely small compared to the statistical uncertainties of 
the measurements. In addition, although the hadronic and partonic calculations yield similar results for 
the correction to the unpolarized cross sections, they yield corrections to the polarization measurements 
with opposite sign at large Q^. The analysis |129j therefore did not to include any TPE corrections to 
the polarization measurement, as either calculation would have had an extremely small impact on the 
final result. 

The extracted form factors from the combined analysis of polarization measurements and TPE- 
corrected cross sections is shown in Fig. [271 Since this publication, there have been updated polarization 
results at high [221 121] and very low [211 [2S1 123 1131] , as well as an extensive set of cross section 
measurements at low [HI]. The global fit of Ref. |129] has been updated in Refs. [231 1131] to include 
the new polarization measurements, and inclusion of the new cross section measurements [81] and a 
detailed evaluation of the uncertainties is in progress |132] . 



5.4.3 Impact on the extracted charge and magnetization radii of the proton 

The proton charge radius is related to the low behavior of the charge form factor. In both electron 
scattering and atomic physics extractions, the root-mean-square radius Re is defined in terms of the 
slope of the form factor at = 0, 



6 



+ 



-6 



(62) 



Because the TPE corrections are finite and have a significant dependence in the limit Q"^ ^ 0, 
as seen in Fig. [131 they can impact the electron scattering extractions of the charge radius. These 
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Figure 27: Proton Gm (left) and Ge (right) form factors from the global fit in Ref. |129] . 
scaled by the dipole form factor fl37|) . Note that above = 6 GeV^ there were no direct 
measurements of G^; and Gm was extracted from the cross section data with an additional 
uncertainty related to the uncertainty in Ge- The solid line is the final fit, while the points 
are from direct extractions from data in small bins. The short-dashed line shows the fit 
from the cross sections with no TPE corrections applied. 



corrections were first included in the second Born approximation by Rosenfelder [50], who found that 
they increase the charge radius by about 0.01 fm. More recently, the full hadronic calculations were 
evaluated |107j and compared to the Coulomb distortion correction [50]. The additional impact of the 
TPE was very small, yielding an additional increase of 0.002 fm. Since Ge dominates the cross section 
at low Q^, it can be extracted reliably without requiring measurements at extremely large scattering 
angles. Therefore, it is not entirely surprising that the charge radius is relatively insensitive to the 
difference between the second Born approximation and the full hadronic TPE calculation which is most 
important at larger scattering angles, as shown in Fig. [T71 

In contrast, the magnetic form factor becomes increasingly difficult to extract from cross section 
measurements at very low Q^. In a Rosenbluth separation, one must extrapolate to e = 0, where the 
reduced cross section is proportional to for — )■ 0. This makes the extraction sensitive to small e- 
dependent corrections, enhancing the sensitivity to the TPE corrections, as well as the difference between 
the second Born approximation and the full calculation. There is a significant difference between two 
recently extracted values of the proton magnetization radius. The first obtains an RMS magnetic 
radius of 0.777±0.017 fm [HI] from a large body of new, high-precision cross section measurements. 
The second finds 0.867±0.020 fm [23] from a fit combining previous cross section measurements with 
new low-Q^ polarization measurements [231 125] . The latter analysis includes hadronic TPE corrections 
for all cross section measurements, while the former includes the Coulomb distortion correction of 
Ref. [3H], corresponding to the = limit of the soft-photon approximation. This is an overestimate 
of the correction at all measured values, and entirely neglects the dependence of the correction, 
which is important in the extraction of the magnetic radius. An estimate of the impact of a more 
complete TPE correction |134] suggests that this could have a large impact on the form factor and 
magnetic radius extracted in Ref. [8T] . 

5.5 Normal asymmetries 

The TPE exchange process gives rise to a nonzero contribution to the elastic cross section for a recoil 
proton polarized normal to the scattering plane. By time reversal invariance this is also equivalent to 
scattering (unpolarized) electrons from a target polarized normal to the scattering plane |133j . Normal 
polarization observables vanishes in the Born approximation, and their measurement directly accesses 
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the imaginary part of the TPE amphtude. A detailed discussion of the formahsm and experimental 
status of normal asymmetry measurements can be found in Ref. ^llOj . 
The target normal asymmetry A at is defined as [llUl 1133^ 1135] 

— , ^ 

An = 7, (63) 

where the cross section for unpolarized electrons scattering from a proton target with spin 

parallel (antiparallel) to the direction normal to the scattering plane, defined by the spin vector Cn = 
k X k' /\kx k'\ (see Eq. f|T2|) )PI At order a in the electromagnetic coupling, the target normal asymmetry 



is given by 



= \Mj- ' ^^'^ 



providing a direct measure of the imaginary part of Ai 
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Figure 28: (Left) Ratio of the TPE normal polarization correction to the unpolarized Born 
contribution as a function of e, for = \ (dotted), 3 (dashed) and 6 GeV^ (solid). (Right) 
Target normal polarization asymmetry, as a percentage, as a function of the center of mass 
scattering angle, 6*011, for Q"^ = 1 (dotted), 3 (dashed) and 6 GeV^ (solid). Figures adapted 
from Ref. [6]- 



In Fig. I2H] (left) the ratio 6n of the TPE normal contribution from nucleon elastic intermediate states 
relative to the unpolarized Born cross section is shown as a function of e for several values of Q^. Note 
that there is no IR contribution to the normal polarization. The normal polarization contribution is very 
small numerically, 6n ^ 1%, and has a very weak e dependence. In contrast to the TPE longitudinal 
and transverse polarization corrections 6l and 6t in Sec. \5.3\ the normal polarization ratio is smallest 
at low e, becoming larger with increasing e. 

The target normal asymmetry in Eq. ( IMj) is shown in Fig. |28] (right) as a function of the center 
of mass scattering angle, O^m, for several values of Q^. The asymmetry is relatively small, of the order 
of 1% at low 9cm for ~ 3 GeV^, but grows with [6]. The magnitude of the elastic contribution to 
An is relatively model independent (see also Refs. [86 |ll33j ). as it is determined mostly by the (on-shell) 
proton electromagnetic form factors, which are reasonably well determined. Estimates of contributions 



^Equivalently, by time reversal invariance, Ajv would be the asymmetry for scattering unpolarized electrons from an 
unpolarized proton target, with the recoil proton polarized normal to the scattering plane. 
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from higher-mass intermediate states have been made within the GPD approach [86], and suggest that 
these may be as significant as the elastic at forward angles ^cm, but very small at backward angles. 

Unfortunately no data exist on target normal asymmetries for elastic scattering on a normal- 
polarized target. Experiments have been performed at Jefferson Lab to measure the target normal 
asymmetry of the neutron on a ^He target for deep-inelastic scattering |136] and quasi-elastic neutron 
knockout [137j . and the data are currently under analysis. Note that the observation of nonzero effects 
here would not necessarily signal TPE effects, as there are likely to be nuclear final state interaction 
effects which can also produce a nonzero asymmetry. A recent search for target normal asymmetries in 
deep-inelastic scattering at HERMES |138j did not find a TPE signal within the ~10^'^ uncertainties 
of the measurement. 

Although not directly relevant to the elastic form factor extraction, the observation of nonzero 
target normal asymmetries would provide direct evidence of TPE in elastic scattering. Knowledge of 
the imaginary part of the TPE amplitude could be used to constrain models of Compton scattering, 
or as input into dispersion relations to obtain the real part of the TPE amplitude from the imaginary 
part (see Sec. 14. 7p . 

The imaginary part of the TPE amplitude can also be accessed by measuring the electron beam 
asymmetry for electrons polarized normal to the scattering plane scattering from unpolarized targets. 
The corresponding beam normal asymmetry is then defined analogously to Eq. ( I63ll . with the 
electron spin parallel or antiparallel to the normal polarization vector (^^r. Since the beam normal 
asymmetry involves flipping the helicity of the electron, is it zero in the limit rrie — )■ 0. Several such 
measurements have been made in elastic ep scattering |139l I140[ 1141] in connection with measurements 
of parity-violating elastic scattering. The most recent result [142j includes both ep and ed scattering 
at backward angles, and the results indicate the importance of including inelastic contributions in the 
intermediate state. 

6 Two-photon exchange in other reactions 

While the effect of TPE has been most dramatically illustrated for extraction of proton electric and 
magnetic form factors using the LT separation method, TPE has also been studied for its effect on the 
extraction of other observables, ranging from elastic neutron and transition form factors, to the pion 
and nuclear form factors. 

6.1 Neutron form factors 

Because the magnitude of the electric form factor of the neutron is relatively small compared to that 
of the proton, especially at low Q^, the effects of TPE may be even more pronounced for than 
for G^^. Within the hadronic framework of Sec. HJ the relative TPE correction to the neutron elastic 
cross section is shown in Fig. [29] as a function of e for = 1, 3 and 6 GeV^ [6], using input neutron 
form factors from Ref. [HS]- For comparison the correction at = 6 GeV^ is also computed with the 
parametrization from Ref. [30], with the difference between these indicative of the model dependence 
of the calculation. 

Since there is no IR-divergent contribution to the TPE correction for the neutron, the total correction 
5" is displayed in Fig. |29] As discussed in Sec. 14.5] for — )■ the TPE correction depends only on 
Fi(0), and is independent of the anomalous magnetic form factor F2. In this limit the TPE correction 
for the neutron therefore vanishes. 

The most notable difference with respect to the proton (Fig. [T3]) is the sign and slope of the TPE 
correction. In particular, the magnitude of the neutron correction is around 3 times smaller than for 
the proton, and the negative slope arises from the negative anomalous magnetic moment of the neutron 
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Figure 29: TPE correction 5" to the unpolarized electron-neutron elastic scattering cross 
section at = 1 (dotted), 3 (dashed) and 6 GeV^ (solid and dot-dashed). The dot-dashed 
curve corresponds to the form factor parametrization of Ref. [30], while the others are from 
Ref. [96]. Figure adapted from Ref. [6]. 



Hn [6|. While the correction for the neutron is smaller than that for the proton, the value of G"^ is small 
at all values, so the effects of TPE on an LT separation will be magnified. Figure [30] (left panel) 
shows the impact that a TPE correction would have on a Rosenbluth extraction of ^nG%/G\j. The 
"uncorrected" values are taken from a global fit [96], and the filled circles show the modification after 
applying the TPE corrections from Fig. [29], assuming that the initial value came from an LT separation 
covering two different e ranges. Inclusion of the TPE effect would improve the significance of the upper 
limits on G^ from the SLAG Rosenbluth extraction |143j up to = 4 GeV^, as well as yielding a small 
shift in GIj. 

While the TPE corrections to the form factor ratio from LT separations are significant at large Q^, 
the neutron G^ form factor is typically extracted from polarization transfer experiments. To compare 
the TPE effects on the ratio fj,nG^/G^ extracted using the polarization transfer method, we take the 
same "uncorrected" starting values for and show the impact of the TPE corrections at 

two values of e, namely e = 0.3 or 0.8 (right panel of Fig. [30]) . The shift in the form factor ratio is 
considerably smaller than that from the LT method, and below the present experimental uncertainties, 
but nevertheless represents a 3 — 4% suppression at = 3 GeV^ and 5 — 10% at = 6 GeV^. 

For kinematics typical of neutron form factor experiments at Jefferson Lab, the ratio G^/G^^ was 
recently measured in experiment E93-038 |144j at = 1.45 GeV^ for e ~ 0.9, at which the TPE 
correction was ^ 2.5%. In the subsequent extension E04-110 |145j at Q'^ ^ 4.3 GeV^, the TPE 
correction for e ~ 0.82 is expected to be around 3%. While small, these corrections will be important 
to take into account in order to achieve precision at the several percent level. Two-photon exchange 
effects will also need to be taken into account when extracting the neutron magnetic form factor G^j 
from cross section data, if precision at the ~ 1 — 2% level is sought. In particular, in measurements of 
the ratio of neutron to proton cross sections, the neutron cross section is obtained by multiplying the 
ratio by the total measured proton cross section which depends on both the Born form factors and TPE 
corrections. The extracted neutron cross section must then be corrected for TPE effects in the neutron 
to obtain the Born form factors. 
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Figure 30: (Left) Effect of TPE on the ratio of neutron form factors fin.G^/G^j using LT 
separation. Tfie uncorrected points (open circles) are sliifted by taking linear fits to the TPE 
correction in Fig. [22] over the quoted e ranges. (Right) TPE effect on polarization transfer 
measurements. The corrected points correspond to e = 0.3 (filled squares) and e = 0.8 (filled 
circles) (offset for clarity). Figure taken from Ref. [6]. 
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6.2 Electroproduction of resonances 

Beyond elastic final states, inelastic electroproduction channels provide additional information on 
hadron structure, from mapping out the spectrum of states in QCD, to probing the spatial distributions 
of hadrons. Of particular interest is the reaction whereby the produced final state is in the A(1232) 
resonance region, which has been recently studied both experimentally [831 1146] and theoretically [147] . 

At the Born level, A electroproduction is parametrized by the same electromagnetic N ^ A tran- 
sition amplitudes that enter the calculation of the A intermediate state contribution to the TPE cor- 
rections in elastic ep scattering in Sec. 14.4.11 It therefore provides an important consistency check on 
the role of the A resonance in electron scattering. Kondratyuk and Blunden |148] computed the TPE 
corrections to the unpolarized cross section for A production in electron-proton collisions, including 
both nucleon and A contributions in the intermediate states in the one-loop diagrams. As in the case 
of elastic scattering, the TPE effects from the intermediate nucleon and A have opposite signs in most 
kinematical regimes, and are pronounced even at low energies. 

It is straightforward to show that the sum of the box and crossed-box amplitudes with an intermedi- 
ate A is gauge invariant by itself. For the nucleon intermediate state, an additional 'j'yNA contact term 
has to be added to the box and crossed-box diagrams. Such a term is required for gauge invariance, 
and can be constructed by the standard procedure of minimal substitution [149] . 

An important theoretical ingredient of the calculation is the 7AA vertex form factor. Various forms 
of this vertex feature prominently in the studies of electromagnetic interactions of the deuteron |150] 
and three- nucleon bound states [151], as well as in the recent extraction |152] of the A magnetic dipole 
moment. 

The TPE corrections in A electroproduction are, as expected, generally much smaller than the Born 
contribution [148] . but can be important for a precise analysis of ep scattering in the resonance region. 
Results for = 1 GeV^ and = 3 GeV^ are shown in Fig. EH Kondratyuk and Blunden [148] find a 
pronounced e dependence of the TPE correction to the value and sign of the (dominant) 7AA magnetic 
coupling constant. A current analysis |83j of experimental Rosenbluth separations in ep scattering 
provides strong constraints on the nonlinearity in e, although additional high-precision data would 
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Figure 31: (Left) Two-photon exchange graphs for the ep — t- eA reaction. The loop diagram 
with the 77A^A contact term (denoted by the black circle) ensures gauge invariance of the 
calculation. (Right) Two-photon exchange correction to the unpolarized cross section for A 
production in electron-proton collisions, calculated for = 1 GeV^ (solid) and = 3 GeV^ 
(dashed). The upper, middle and lower pairs of lines are labeled by the values of the 7AA 
coupling constant ("N" corresponding to = 0, i.e. to the absence of an intermediate 
A state). Figures adapted from Ref. |148] . 



allow more definite conclusions to be reached. 

For further progress in the evaluation of higher-order effects in electron-nucleon collisions, a more 
detailed knowledge of the 7AA vertex is needed. Being formulated in terms of hadronic degrees of 
freedom, this model is somewhat complementary to the approach of Ref. jl47j . where TPE effects were 
calculated using the formalism of generalized parton distributions. 

6.3 Timelike form factors 

Electromagnetic form factors in the space-like {Q^ > 0) and time-like (Q^ < 0) regions both yield 
information on the structure of hadrons. The form factors are real in the space-like region, while 
the time-like form factors have a phase structure reflecting the flnal-state interactions of the outgoing 
hadrons, and are therefore complex. 

The interference of one- and two-photon exchange can be studied in the process e~^e~ — ?■ pp, and 
is related by crossing symmetry to the elastic ep interaction. The interference manifests as an angular 
asymmetry in the differential cross section. This was studied recently by Chen et al. |153j in a hadronic 
model with essentially the same ingredients as the ep calculations. The corrections are again at the 
few percent level, and are shown in Fig. [32] for = 4 GeV^. A search for TPE effects in BABAR 
e'^e~ — )■ pp data |154j showed no indication of the expected forward-backward asymmetry within the 
precision of the data |155] . 

In addition to the unpolarized differential cross section, Chen et al. |153] considered the single-spin 
polarization observable Py, and the double-spin polarization observables Px and P2, which require the 
incoming electron to be polarized. They suggested that Pz in particular should be considered in future 
experiments looking for TPE effects. 

6.4 Pion form factor 

As the lightest bound state of quarks and antiquarks, the pion plays a unique role in QCD, and 
determining its electromagnetic structure is of great importance for understanding the realization of 
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chiral symmetry in nature. Extractions of the pion form factor in the space-hke region from 

measurements of the pion electroproduction reaction ep — >■ evr+n have recently provided data on the 
Q"^ dependence of up to values of ~ 2.5 GeV^ |156[ I157[ 1158] . and higher measurements are 
planned to ~ 6 GeV^. The role of TPE in electromagnetic scattering from the pion was recently 
investigated by several authors P 11171 11591 1160] . The analysis of TPE from a spin-0 target is, in fact, 
considerably simpler than that for spin-1/2 targets. 

The form factor of the pion is defined through the matrix element of the pion current, 

{7,{p')\J^m7r{p)) = {p + pYF^iQ'). (65) 

In the Born approximation the amplitude for an electron scattering from a Tr"*" is given by 

M; = -j^u,{k')^^u,{k){p + p'YFM), (66) 

and the Born cross section in the target rest frame is 



cTMottF^iQ'), (67) 



where dMott is given in Eq. (p. 

The TPE amplitude for the box diagram has the form [HI 1117] 



= -te'J ^ ViJrAF(gi,A)A^(g2,A), (6J 



where the leptonic tensor Lf^,y is given in Eq. (!32|) . and the crossed-box is obtained using crossing sym- 
metry. The pion hadronic tensor i?^'^ in principle contains contributions from all hadronic excitations 
in the intermediate state. For the dominant pion elastic contribution, one has 

= {2p + qi + qrF^{Ql)AF{p + q,,m^){2p + q^yF^{Ql). (69) 

The pion form factor is then modified according to F^{Q^) — )■ F^{Q^) (1 + 6'^), where, as for the 
nucleon TPE correction in Eq. (l35|l . the relative correction 6'^ is given by 
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Figure 33: Two-photon exchange correction 5 to the pion form factor squared, relative to 
the Mo-Tsai (MoT) contribution [3Z1 EB], as a function of e for various (left), and of 
for various e (right). Figures adapted from Ref. [9j. 



As in the case of the nucleon electromagnetic form factors, experimental analyses of pion form factor 
data typically use radiative corrections computed by Mo and Tsai in the soft-photon approximation 
[33, ES] (see Sec. 14.11) . The effect of the IR-finite, structure-dependent contribution is illustrated in 
Fig. [331 where the difference between the full TPE correction and the Mo-Tsai prescription [37t [38] , 
5 = 8"^ — (5ir(MoT), is shown as a function of e and Q^. Here a monopole parametrization for the 
"bare" pion form factor in Eq. fl68|) is used, 

F.m = (i+f^) \ (71) 

with the cut-off parameter A^^ = 770 MeV corresponding to the p-meson mass. The sensitivity of the 
results to the input pion form factor was studied in Ref. [9j, using an alternative parametrization |161] 
which gives a better description of the available data and builds in gauge invariance constraints for the 

— )■ limit and perturbative QCD expectations for the — t- oo behavior. The differences at low 
are negligible, but become noticeable at high Q^, although do not affect the results qualitatively. 

At low (Q^ ~ 0.01 GeV^) the TPE correction is positive and of the order of 1% at backward 
angles (small e), decreasing to zero as e — )■ 1. For larger the correction becomes more negative up 
to ~ 1 ~ 2 GeV^, especially at backward angles, but changes sign at intermediate e. Unlike for ep 
scattering, the en cross section vanishes in the extreme backward limit (e: — ?■ 0). Above ~ 2 GeV^ 
the correction grows once again, reaching ~ 1% at = 10 GeV^ [H]. Interestingly, the TPE correction 
is most positive at very small {Q"^ <^ 1 GeV^) and at large {Q"^ ^ 1 GeV^), reaching its minimum 
values at ^ 1 - 2 GeV^. 

The contributions from higher-mass intermediate states to the pion hadronic tensor were considered 
in Ref. |117] within a dispersion relations approach. Because the mass difference between the pion 
and the next excited resonant state, the p meson, is ~ 5 times as large as the pion mass, one would 
not expect large contributions from excited hadronic intermediate states. This was indeed confirmed 
in explicit calculations of the p and 6i(1235) meson contributions, which were found to be negligible 
for < 2 GeV^, as Fig. |3l] illustrates. At larger (> 4 GeV^) the inelastic contributions become 
comparable to the elastic for e ^ 0.2; however, in the region at low e where the TPE effect is greatest, 
these are still significantly smaller than the elastic components. 

In contrast to the proton form factor case, where the TPE effects give large corrections to the elastic 
form factors extracted from LT separated cross sections at large 1129], the TPE corrections to the 
pion form factor are relatively small. This stems from the fact that electron scattering from a scalar 
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Figure 34: Elastic and inelastic (p and 6i(1235)) contributions to the TPE correction 6 
(in percent) relative to the Maximon-Tjon result ^49j, as a function of e for = 1, 2 and 
4 GeV^. Figure adapted from Ref. |117] . 



target is described by a single form factor, with no LT separation necessary. On the other hand, since 
is extracted via LT separation of the pion electroproduction cross section, TPE with one photon 
attached to the pion and the other to the initial proton or final neutron could modify the longitudinal 
cross section, and may need to be considered. Furthermore, the effects on the pion form factor from 
nonresonant contributions in the intermediate state have not yet been evaluated, and may also need to 
be considered in future analyses. 



6.5 Electron-nucleus scattering 
6.5.1 TPE in deuteron form factors 

The TPE corrections to the (spin-1) deuteron (D) elastic electromagnetic form factors and to the 
e~^e^ —7- DD process have been discussed in Refs. |162l I163[ 11641 1165[ 1166] . The reaction amplitude 
contains six generalized form factors, but only three linearly independent combinations of them (gener- 
alized charge G^, quadrupole Gq, and magnetic form factors) contribute to the cross section |162] . 

Dong and collaborators |163[ I164[ 1166] estimated all the TPE corrections to the conventional form 
factors of the deuteron using an effective Lagrangian approach. In this model, the two photons couple 
to one of the two nucleons in the deuteron. Dong found the TPE corrections to be small (less than 
1% for < 2 GeV^), and most significant for G^ |163] . He also suggested that two of the additional 
form factors could be tested in measurements of the double and single polarization observables |163] . 
The effect makes a sizable contribution at backward angles to the polarization observable Py of eD 
scattering, which vanishes in the Born approximation |166] . 

Kobushkin et al. fl6^ also discuss elastic electron-deuteron scattering beyond the Born approxima- 
tion, including contributions where the two photons can also interact with different nucleons. They 
conclude that TPE may give a large contribution to elastic eD scattering, but point out that their 
estimates have large uncertainties. The most important source of uncertainty comes from the short- 
range part of the deuteron wave function. It is suggested that experimental study of TPE in elastic eD 
scattering for in the few GeV^ range can give important information about the deuteron structure 
at short distances. 
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6.5.2 ^He elastic form factors 

The hadronic formalism of Sec. Hllias also been applied to the case of elastic scattering from the spin-1/2 
^He nucleus [6], where contributions of unexcited ^He intermediate states were computed. Of course 
these are likely to constitute only a part of the entire TPE effect, as contributions from break-up 
channels may also be important. Nevertheless, a comparison of the size and magnitude of the TPE 
corrections from elastic intermediate states in ^He relative to those in the proton illustrates a number 
of features of TPE. 

The expressions used to evaluate the TPE contributions to the ^He correction S^^'^ are similar to 
those for the nucleon in Sec. 14.31 but with some important differences. The charge of the nucleus is 
Ze, with Z = 2 for ^He, and the mass MsHe is approximately 3 times larger than the nucleon mass, 
while the anomalous magnetic moment is KSHe = —4.185. The internal photon-^He form factors are also 
softer than the corresponding proton form factor (the charge radius of the ^He nucleus is ~ 1.88 fm), 
and have zeros at ~ 0.45 and 0.7 GeV^ for the charge and magnetic form factors, respectively jl67] . 



0.04 



u" 

-0.04 



-0.08 - 



-0.12 




Figure 35: TPE contribution to the unpolarized electron-'^He cross section, with the ^He 
elastic intermediate state, as a function of e, for = 0.05 (sohd), 0.2 (dot-dashed), 0.5 
(dashed) and 1 GeV^ (dotted). A parametrization of the form factor from Ref. [167j is used 
in all cases, except for the upper solid and dot-dashed curves, which uses a dipole with mass 
AsHe = 0.37 GeV. Figure adapted from Ref. [B]. 



The TPE correction S , relative to the Mo-Tsai IR-divergent result, is shown in Fig. |35] as a 
function of e for several values of Q^, using the ^He form factors from Ref. [167J. As for the proton, 
the correction is negative at low Q^, and generally increases in magnitude with increasing Q^. The 

dramatic change of sign in 6 at Q"^ = 0.5 GeV^ stems from the presence of the zeros in the form 
factors at ~ 0.5 GeV^. 

-3 He 

To estimate the model dependence of the results, Blunden et al. also calculated 5 assuming a 
dipole shape with a cut-off mass AsHe = 0.37 GeV, fitted to the ^He radius, which gives a reasonable 
approximation at low {<, 0.05 GeV^). The results with this form factor are about a factor 2 
smaller in magnitude than for the form factor from Ref. |167] . At larger Q"^ the dipole shape is a less 
reliable representation of the ^He form factor, making it difficult to estimate the model dependence of 
5 The results in Fig. |35] illustrate the potential relevance of TPE effects for future ^He form factor 
measurements, which may reveal interesting TPE effects at larger |168j . 
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7 Parity-violating electron scattering 



In addition to the exchange of one or more virtual photons between the electron and nucleon, the 
Standard Model allows the scattering to take place via the exchange of a neutral Z boson. Because the 
Z boson mass is some two orders of magnitude larger than the proton mass, the weak exchange process 
is strongly suppressed relative to the electromagnetic reaction. Nevertheless, asymmetries sensitive to 
the 7Z interference amplitude, which are of the order of several parts per million (ppm), have been 
measured in modern accelerator facilities. These have been used to probe the strangeness content of 
the proton, through measurements of the strange electric and magnetic form factors |169[ 11701 11711 1172[ 
11731 11741 11751 11761 11771 1178j , as well as to measure the weak charge of the proton [179j . 

The ^Z interference term is isolated by polarizing the incident electron and measuring the difference 
between right- and left-handed electrons scattering from unpolarized protons. A parity-violating (PV) 
asymmetry can then be defined in terms of the differential cross sections as 

^pv = ^ ^ , (72) 

where is the cross section for a right-hand (helicity /i = +1) or left-hand (helicity /i = — 1) polarized 
electron. The numerator in the asymmetry is sensitive to the interference of the vector and axial- vector 
currents, and hence violates parity. The natural size of the PV asymmetry is ~ Q^/Mf, which for 

~ 0.1 GeV^ is of the order 10~^ — 10~^ — values that can be routinely measured at current facilities 
such as Jefferson Lab or MAMI at Mainz |180j . 

In view of the large TPE effects on electromagnetic form factors discussed in the proceeding sections, 
the question naturally arises of what effect the exchange of two bosons (7 or Z) may have on PV 
asymmetries. Because both the strange form factors and the proton weak charge are numerically small 
quantities, these contributions could affect their extraction significantly. 

One should note that the "standard electroweak radiative corrections" already include a contribution 
from two-boson exchange (TBE) effects. These are usually taken from the classic calculation of Marciano 
and Sirlin |18H 1182] . and include 7Z, ZZ, and WW exchange box and crossed-box diagrams. Because 
of the large masses of the Z and W bosons, the box diagrams contain contributions for virtual four- 
momenta of all mass scales. For the ZZ and WW boxes, the high mass scales dominate, and a 
calculation in terms of the quark structure of the nucleon is reliable. However, for the 7Z diagram both 
high and low mass scales contribute. The high-mass scale, calculated in terms of quarks, accounts for 
most of the effect, but the low momentum contribution in terms of a hadronic picture is not negligible. 
The quark-based calculation is discussed in more detail in Sec. 17.31 

In this section we review the TBE corrections to PV asymmetries arising from the interference 
between single Z boson and 77 exchange amplitudes (which we denote by "^(77)"), and between 
the one-photon exchange and 7Z interference amplitudes (denoted by "7(^7)"). Because PV electron 
scattering experiments are typically performed at values ^ 1 GeV^, the hadronic formalism described 
in Sec. HJis the more natural implementation of hadronic structure effects. For a detailed examination 
of the GPD-based approach see Refs. |110[ 1135] . 



7.1 Parity-violating asymmetries 

In the Born approximation, the weak neutral current amplitude is given by 

where g = e/ sin^vK is the weak coupling constant, Mz is the Z boson mass, and the strength of the 
effective four-fermion interaction is given by the Fermi constant Gp = 7ra/(-\/2M| sin^ Ow cos^ Ow)- At 
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tree level the weak mixing angle 9w is related to the weak boson masses by siia^ 9w = 1 — M^/M^, 
where Mw is the W boson mass. The weak leptonic current is given by a sum of vector and axial-vector 
terms |183] . 

jz^ = Ue{k',h)-ff,{gy - g'X-f5)ue{k,h) , (74) 
where the vector (gy) and axial- vector {g'^) couplings of the electron to the weak current are defined as 

9v = -\{l-4smHw) , 9l = -\, (75) 

respectively. Note that some of the TPE hterature P [951 UnOl [1801 [M HHHl [1861 [IHH HHHl HHH] uses 
definitions whereby the vector and axial charges are scaled by a factor 2 relative to those in Eq. ( j75l) . 
so that care should be taken when comparing formulas for amplitudes. 
The matrix elements of the weak hadronic neutral current is given by 

Jz = UN{p')T''z{q)uN{p) , (76) 
where the current operator is parametrized by three weak nucleon form factors, 

r^(?) = YF^iQ') + '-^FriQ') + Yi^criQ') . (77) 

Alternatively, the current operator can be expressed in terms of the Sachs weak form factors Gf^/, in 
analogy with the electromagnetic nucleon form factors in Eq. ( ITTj) . 

GT(g2) = F^iQ') - rFi^m , Gl^{Q') = F.'^Q^) + F^^m . (78) 

The differential cross section including 7 and Z exchange is given by Eq. (|8]), where the amplitude 
is now a sum of the 7 and Z Born contributions, 

\M\^ = \M^ + Mzf = \M^\^ + 2ne{M;Mz) + \Mz\^ ■ (79) 

At the kinematics of interest in this review, the purely weak contribution |A^^p is small compared with 
the other terms and can be neglected. The purely electromagnetic contribution cancels in the numerator 
of the asymmetry in Eq. ( 1721) . so that Apy is sensitive to the parity- violating part of 27le (A^*A^^), 
involving the interference of with the product of vector and axial-vector currents in A4z- The 
vector- vector and axial- axial parts of A^^ also cancel in the asymmetry. The denominator is dominated 
by the electromagnetic term, |A^-yp. Evaluating the asymmetry explicitly, one has for PV asymmetry 
of the nucleon 

<v = - ^ {-29l {sGfG^ + r<Gff ) + 2,^.'<Gr } , (80) 



where e' = a/ t(1 + r){l — e'^), and cr]^'^ is the reduced 7A^ cross section defined in Eq. (jTOl) . 

For a proton target, assuming isospin symmetry (see, however, Refs. |190i I191j ). the weak vector 
form factors G^^^^ can be related to the electromagnetic form factors of the proton (neutron) G^^^"^ by 

Ge^m = (1-4 sin^ 9w) G^^^ - G^^^ - G% j^^, (81) 
where G% ^.j are the contributions from strange quarks, and 

G^^(O) = = l-Asin^Ow [Born approximation] (82) 
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is the weak charge of the proton. Because Q^y is numerically small, the overall contribution to G^j^ 
from the proton electromagnetic form factors is suppressed. The weak axial-vector form factor of the 
proton is given by = —G^ + G\, where 6*^(0) = = 1.267 is the axial-vector charge, and G\ is 
the axial-vector strange quark contribution, related to the spin of the proton carried by strange quarks 
as measured in deep-inelastic scattering |192[ 1193] . In terms of the proton vector {Ay), axial- vector 
[Aa) and strange (As) contributions, the asymmetry can then be written as 



where 
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where is the reduced cross section for a proton target. 



7.2 Two-boson exchange corrections 

Beyond the Born approximation, the PV asymmetry, Eq. ( IHOjl . receives corrections from higher-order 
radiative effects, including both the electromagnetic two-photon exchange contributions discussed in 
Sec. m and corrections involving 7-Z boson loops. There are several ways in which the PV asymmetry 
can be represented in the presence of higher-order radiative corrections. The approach pioneered by 
Marciano and Sirlin |181l 1182] parametrizes the electroweak radiative effects in terms of parameters p 
and K, such that the weak charge of the proton in the presence of higher-order corrections becomes 



w 



p (l - An sin^ 9w] 



i7) 



Including the standard radiative corrections, such as vacuum polarization, vertex corrections, and 
hadron structure-independent two-boson exchange contributions, one has p = 0.9877 and k = 1.0026 
|183] . The TBE amplitudes discussed here give additional corrections to p and k, which are denoted in 
Refs. Unni MS\ M by Ap and Ak, respectively. 

An alternative parametrization is in terms of isoscalar and isovector weak radiative corrections for 
the vector form factors, and a similar set of corrections for the axial-vector form factors. In this case 
the vector part of the PV asymmetry is written 



A 



V 



1 
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(88) 



where the proton and neutron radiative corrections are given, to first order in p — 1 and k — 1, by 
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The radiative corrections to the strange part of the asymmetry (|86l) enter via a multiplicative factor 



A, Asil + R 



(0)^ 



(90) 
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Figure 36: An example of a two-boson exchange contribution to parity-violating elastic 
scattering, due to 7Z exchange (given by the amplitude A^^z)- Other contributions include 
the amplitude with the 7 and Z interchanged (denoted by A^^^), and the crossed-box 
diagrams. 



where the isoscalar radiative correction i?y^ = p — 1. For the axial asymmetry A a-, the generalization 
of the tree level axial form factor, — )■ , implicitly contains higher-order radiative corrections for 
the proton axial current, as well as the hadronic anapole contributions [1841 1194] . 

More explicitly, the corrections to the cross section arising from the 77 and 7Z TBE contributions 
can be obtained from Eq. (1791) by the replacements 

M^ + M^^, Mz ^ Mz + M^z + Mz-y , (91) 

where Ai^z and Aiz-y correspond to the two different orderings of the 7 and Z exchanges in the box 
and crossed-box diagrams. Fig. [361 

Because the coupling of the Z boson to the electron and to the proton contains both vector and axial- 
vector components, the interference amplitude M^z (and similarly Aiz-y) contains (parity-conserving) 
vector electron-vector hadron and axial electron-axial hadron terms from the Z exchange, along with 
(parity-violating) vector electron-axial hadron and axial electron-vector hadron terms (the 7 exchange 
contributions are pure of course pure vector- vector). The parity- violating asymmetry Apy selects only 
the vector-axial vector part of the amplitude, Ai^^^\ which can be decomposed into the vector hadron 

contribution A4^^^^ and the axial- vector contribution Ai^^^K These have the crossing symmetry 
properties under the interchange s ^ u {ot E -H- —E'), 

M'^r^'^-M = -M';r'-{s,t) , (92) 
M^r^'-'^^'M = +M^r'-{s,t) , (93) 

and similarly for 

The relative corrections from the 2^(77) and 7(7^^) interference terms to the PV asymmetry can 
then be identified as 

5zi,,) = f , (94) 

27^e (mIM^^'P + MlMfP) 
S.iz.) = ^ 7^ ,pv)X ^ > (95) 



with the full PV asymmetry, including TBE corrections, given by 



= (1 + S)A% ^ { + ] , (96) 
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where Apy is the Bom level asymmetry in Eq. (|80|) . The electromagnetic TPE correction 57(77) froi^ 
Eq. f l5B]) is typically only a few percent [SI El E] , so that the full correction 6 can be written approximately 
as 

S ~ '^^{77) + "^7(^7) ~ "^7(77) • (9''') 

The amplitudes Ai^^, M.-yz and M-z-y contain contributions from both nucleon elastic and inelastic 
intermediate states. At forward scattering angles the contributions from the high-mass inelastic contin- 
uum can be described with the help of dispersion relations, discussed in Sec. l7.3l below. For non- forward 
angles, the elastic and inelastic terms must be computed directly, along the lines of the 77 corrections 
in Sec. ID 



7.2.1 Nucleon elastic contributions 

The calculation of the 7Z correction with a nucleon intermediate state to the ep ^ ep cross section, 
where the incident electron is polarized with helicity h, mirrors that of the 77 amplitude Ai^-y in 
Eq. dnni), with the replacement of the -fNN vertex function by in Eq. ([77]) HDHl IISS], 

M^^z = / ^ Km ^^(^1' ^) ^^(^2, Mz). (98) 

Here the leptonic tensor L^f for the box diagram is given by 

Liu = Ueik', h) {gy-f^ - gA7f^75) Spik - gi, me) 7^ Ueik, h). (99) 
The hadronic 'jZ interference tensor Hj^^^-^j^ for the nucleon elastic intermediate state is given by 

K^z)N = UN{p)^'ziQ2)SF{p + q,,M)T;iq^)uNip). (100) 

Similar expressions hold for the conjugate amplitude J^z-y Only the PV parts of the amplitudes Aiz'y 
and Ai^z contribute to the PV asymmetry Apy. 

While the electromagnetic nucleon form factors used in the jNN vertex functions are relatively 
well constrained by data, for the weak ZNN vertex the form factors are generally less well determined. 
For the weak vector current, the form factors can be directly related to the 'yNN form factors through 
the conservation of the vector current (CVC). The axial- vector current, on the other hand, is not 
conserved, however, some constraints on the axial form factor have been obtained from elastic neutrino 
scattering data. Using a dipole fit, G^^(Q^) = —(7^/(1 + (5^/A^(-^-))^, one finds that the mass parameter 
^Ar(yi) ~ ^ GeV. Since one of the main purposes of the PV experiments is to extract strange contributions 
to form factors by comparing the measured asymmetry with the predicted zero-strangeness asymmetry, 
in the computation of the TBE corrections one can set the strange form factors to zero, = = G^. 

The contributions from the nucleon intermediate states to the TBE correction 6^ are illustrated in 
Fig. [37] as a function of e for = 0.01 and 1 GeV^. At low the 7(77) and ^'(77) contributions 
are very similar, and therefore partially cancel in the asymmetry Apy. Consequently the asymmetry is 
determined mostly by the 7(2'7) component. At larger (> 1 GeV^) the 7(2'7) component decreases 
in magnitude, while the 7(77) and ^'(77) pieces become large and more negative O [6l |95]. The 
dependence of the total correction 5^ on the input form factors was found P, to be very weak for 
all e, and only becomes appreciable at large (Q^ ^ 1 GeV^). The correction at Q"^ = 0.01 GeV^ 
also appears relatively flat over the range 0.1 <, e ^ 0.8, before dropping rapidly as e — )■ 1. At large 

the total TBE correction becomes more strongly e dependent, decreasing in magnitude at forward 
scattering angles but increasing at backward angles {e — )■ 0). 
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Figure 37: TBE corrections with nucleon intermediate states, for the 7(77) (dotted), 2' (77) 
(dashed) and 7(^7) (sohd) contributions at = 0.01 (left) and 1 GeV^ (right). The 

—N 

finite part of the correction is defined with respect to the Mo-Tsai IR contribution, 6 = 
^^~^ir(MoT) [371 |38]. The 7(77) correction enters with the opposite sign in the asymmetry, 
Eq. ( !97l) . Figure adapted from Ref. [95] . 



7.2.2 A intermediate states 

Contributions to TBE amphtudes from the excitation of the A(1232) resonance were recently considered 
by Tjon et al. j95j and Nagata et al. |186] . Evaluation of the relevant TBE amplitudes follows that of 
the electromagnetic contribution in Sec. 14.4.11 extending the formalism of Kondratyuk et al. [7j to the 
weak sector. 

For the ZNA vertex both vector and axial-vector contributions enter. The vector transitions are 
required by CVC to have the same form as the 7A^A in Eq. ( l38l) . 

TzT^AiPA^q) = ^^{9r{Q')[9''WA-^iyi-ri'q-PA+^Ai'q"] 

+ qUQ") [O^ - g'^'Q ■ pa] + [q' (7>^ - ^?"Va) + mA - rq ■ pa)] ]i,m) 

where again the factor y^2/3 is associated with the — )■ A weak isospin transition. Using CVC and 
isospin symmetry, the vector ZNA form factors can be related to the 7A^A form factors by |95] 

gY{,Q^) = 2{l-2^inHw)9i{Q^). ^ = 1,2,3 (102) 

where the dependence of the electromagnetic 'jNA form factor is parametrized as in Sec. I4.4.1I 

For the axial-vector vertex, nonconservation of the axial-vector current implies the existence of an 
addition form factor. One can use the partially conserved axial current (PCAC) hypothesis to relate 
two of the form factors, leaving a similar expression to that in Eq. f llOip . 

rrilA(PA,g) = ^^{9t{Q')[9''WA-^iy'k-i"rq-PA+^AYq"] 

+ g^iQ') - 9"'q ■ pa] + [q' {TpI - ^?">a) + q' (^Va - Tq ■ pa)] ]■ (i03) 

Note that here the weak isospin transition factor has been absorbed into the definition of the couplings 
[1951 [T96] . 
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Figure 38: Total finite parts of tlie TBE corrections 5, relative to the Mo-Tsai contribution 
[371 EH] 5 with nucleon (dashed) and A (dotted) intermediate states, as well as the sum (solid), 
at = 0.01 GeV^ (left) and 0.1 GeV^ (right). Figure adapted from Ref. [95]. 



The axial form factors are less well determined, but some constraints have been extracted from 
analysis of v scattering data. In a recent analysis, Lalakulich et al. |195[ 1196] parametrized the 
vN — )■ /iA cross sections from bubble chamber experiments at low in terms of phenomenological 
form factors. The available data can be described by the form factors |95l [1961 ttH] gt{Q'^) = 0, 
g^{Q^) = (Q'^/AM^) g^{Q^), with the couplings at = determined in Ref. [I96], and the 
dependence given by a dipole form having a cut-off mass of Aa(a) = 1-0 GeV. Note that Ref. |186] uses 
a different basis of form factors than in Eqs. fllOip and f ll03p : for the relation between these see [93] . 
As for the electromagnetic case, the vertex with an incoming A can be obtained from Eq. ( 139|) . 

The interference amplitude for the box diagram with a A intermediate state can then be written 
[951 [T86] 

M^, = -ze' LlfH^;,^^AF{quO)Aj.{q,,Mz), (104) 

where the leptonic tensor L^^ is the same as in Eq. f l99|) . The hadronic interference tensor Hj^^^-^^ for 
the A intermediate state is now 

^(7Z)A = UN{p')T''^^z^{p + qi,-q2) Sai3{p + qi,M^)T^''^^^{p + qi,qi)uN{p), (105) 

where F^'^^^^ is the sum of the vector and axial- vector vertices in Eqs. fllOip and f llOSp . The corre- 
sponding amplitude Ai^z can be derived in a similar manner. 

The total TBE contribution from A intermediate states, relative to the Mo-Tsai IR result [371 [38] . 
is shown in Fig. [381 as a function of e for = 0.01 and 0.1 GeV^. At low values the A contribution 
is dominated by the jiZj) term. The negligible 77 interference with the Born 7 or Z exchange reflects 
the vanishing of the two-photon exchange correction in the — ?■ limit. For comparison we also show 
in Fig. [38lthe nucleon elastic contribution, and the sum of the nucleon and A terms. The A correction 
is strongly suppressed at low e, but grows with increasing e, becoming as important as the nucleon 
elastic part near the forward limit, e — )■ 1. It is interesting to contrast this with the role of the A in 
TPE corrections to electromagnetic scattering discussed in Sec. 14.4.11 where the A contribution was 
negligible at forward angles but becomes important at backward angles, where magnetic scattering is 
dominant. 

The increased magnitude of the A contribution to PVES at forward angles reflects the growth of the 
invariant center of mass energy for fixed as e — ?• 1 [951 1186] — since the A intermediate state am- 
plitudes A^fty and Ai^z have numerators with higher powers of loop momenta than the corresponding 
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nucleon amplitudes, 6^ will grow faster with invariant energy than 5^. One method to ensure consis- 
tency with the physical bounds at asymptotically high energy requires is to use so-called "sideways" 
form factors, which depend on the center of mass energy as well as on the mass of the virtual boson. 
Alternatively, a dispersive approach, where the intermediate states are on-shell, may be adopted, such 
as that at forward angles discussed in Sec. 17.31 In the absence of a detailed analysis of the high-energy 
behavior of the higher-mass contributions, the predictions of the A calculation should not be taken too 
literally at very large e ~ 1. 

7.2.3 Effects on the strange form factors 

While the TBE contributions to the parity-violating asymmetry in elastic ep scattering are relatively 
small, there are two ways in which the impact on the measurements can be enhanced. Although the 
asymmetry does not have large corrections, these measurements are typically looking for contributions 
from strange quarks which are seen through very small deviations from the asymmetry expected in 
the absence of strange quarks. A few-percent correction to Apy will therefore be significant if the 
contribution from strange quarks is at the level of 5 — 10%. In addition, the corrections discussed in 
the previous section show the difference between the Born prediction for Apy and the value including 
TBE contributions, assuming that one is starting from the correct electromagnetic form factors. There 
is an additional effect for measurements which used proton form factors extracted without accounting 
for TPE corrections in determining the value of Apy for the case with no strange quark contributions. 

In Ref. |13U] . the impact of TPE on the electromagnetic form factors and the TPE corrections 
to the parity-violating asymmetry were examined, but 'jZ box contributions were not included. The 
direct TPE contributions to the asymmetry were found to be ^ 1% for < 3 GeV^, and smallest 
at large e, where the highest precision measurements were planned. However, the impact of using 
proton form factors that did not have TPE corrections applied was found to be much larger, well 
over 5% at 3 GeV^, with the largest corrections in high-e region where high precision measurements 
were planned. A prescription was given that allowed for approximate corrections to be applied, without 
explicit inclusion of the TPE corrections. But because this estimate did not include the 7Z contributions 
and the TPE-corrected form factors can now be reliably extracted, this procedure is now most relevant if 
looking at previous extractions of the strangeness contribution based on TPE-uncorrected form factors. 
However, the use of TPE-uncorrected form factors yields the largest TPE-related error in most of the 
earlier extractions, as the absolute corrections for the two-photon and 7Z diagrams is relatively small. 
However, because these corrections modify all of the extracted strange form factor contributions in a 
highly-correlated fashion, even these small corrections may become important in global analyses of the 
strange quark contributions. 

The full effects of the TBE corrections (two-photon and 7Z exchange) on the Apy asymmetry at 
kinematics corresponding to experiments [IMl [HB ESI [l73l MM MB MM EH EH] designed to 
measure the strange quark form factors of the nucleon were considered in Refs. [951 HOOl I185[ 1186] . For 
the forward angle HAPPEX |169j and GO |173j measurements, the nucleon correction 6^ was found to 
be in the vicinity of ~ 0.1 — 0.2%, increasing to ~ 1.0 — 1.5% for the backward angle GO |174] and 
the earlier SAMPLE \175\ measurements. In contrast, the A contribution 6^ is almost negligible at 
backward angles, but becomes more important at forward angles, where as discussed above its evaluation 
is more questionable. 

The impact of these corrections on the strange form factors is difficult to gauge without performing 
a full reanalysis of the data, since in general different electroweak parameters and form factors are 
used in the various experiments. Nevertheless, following Zhou et al. |185] attempts have been made 
[H I100[ 1186] to quantify the effect from the TBE corrections on the combination G% + PG%j measured 
in the PV experiments, where /3 depends on the kinematics of a particular experiment. Zhou et al. 
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assume the experimental asymmetry A^y can be written as |185] 

A'^^ip, k) = Aly{p\ k'){1 + 5^ + 5^), (106) 

where Apy is the Born asymmetry flHOj) . and p' = p — Ap and k' = k — An remove the existing 
(hadron structure-independent) two-boson exchange corrections contained in p and k [181\ 1182] . From 
the calculated TBE corrections one can then determine Apy and extract the strange asymmetry Ag fl5^ . 
Defining the relative correction 6g to the strange form factors by G^+f3Gf^,^ = {G'^^+ f3G\.f) {1+6g) , where 
G^ is extracted from the experimental asymmetry and G|; is obtained from the Born asymmetry, 
one has |185j 

1 r ,exp AAp \ , GfQ^ ■ 2n . . 



< Ap^ d H 1= — psm OwAk - Aa — } ■ (107) 



Because the magnitude of the strange form factors is generally rather small |194[ 1198] , dividing by the 
strange asymmetry in Eq. (11071) can give a quite large relative correction Sg- The size of the correction 
was found to be typically within 10-20% for most PV experiments |100] . 

The model dependence of these corrections was explored in Ref . [8] , where an estimate was made of 
the induced difference in the strange asymmetry extracted using various input electroweak form factors. 
Comparing results with empirical |129] and monopole form factors (as used in Ref. |185] ). the effects 
ranged from around 15% for the HAPPEX kinematics |169[ 1172] to over 30% for the PVA4 kinematics 
|176] . One should caution, however, that these values are indicative only, and a more detailed reanalysis 
of the strange form factor data including TBE effects a priori would be needed to reach more robust 
conclusions. 

In addition, the fractional contributions may not reflect the importance of the corrections because 
the strange quark contributions are related to the difference between the measured asymmetry and the 
expected asymmetry in the absence of strange quark contributions. Thus, a large fractional correction 
could simply indicate that the strangeness contribution is extremely small. It may therefore be more 
useful to compare the size of the TBE corrections to the uncertainty in the Apy measurements. 

A preliminary attempt at incorporating TBE corrections in a global strange form factors analysis 
was made recently by Young [199] . Using PV scattering data below = 0.3 GeV^ the TBE corrections 
in the proton were found to modify the strange electric and magnetic form factors by 

= 0.0025(182) — > 0.0023(182), = -0.011(254) — > -0.020(254), (108) 

at a scale = 0.1 GeV^. The effect on the strange magnetic form factor is an almost factor two increase 
in the magnitude; while signiflcant, this is still well within the current experimental uncertainty. The 
shift in the strange electric form factor is somewhat smaller. This is mostly because the G^ form factor 
is determined largely by the ^He data, and the TBE effects here have not yet been computed. Overall, 
though, the conclusion appears to be that TBE effects provide relatively mild corrections to strange 
quark form factors. On the other hand, a signiflcantly larger effect from TBE has been found in nearly 
forward PV electron scattering at very low Q^, which we discuss next. 



7.3 7Z corrections to the proton weak charge 



In parallel with the discussion of Sec. 14.71 we consider the forward scattering limit, where k' = k — q, 
with t = kept small but flnite. The parity- violating proton asymmetry in Eq. ([HSD is related to the 
weak charge of the proton Qy^ [184] . 
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Including electroweak radiative corrections, the proton weak charge is defined at zero electron energy 
E and zero momentum transfer as [200j 

= (1 + Ap + Ae)(l-4sin2^^(0) + A^) + a^w + ^zz + D-yzW, (HO) 

where sin^6'vi/(0) = 0.23867(16) is the weak mixing angle at zero momentum, and the corrections Ap, 
Ae and Ag are given in Ref. |200j . The contributions ^ww and D^z from the WW and ZZ box and 
crossed-box diagrams can be computed perturbatively, while the 7Z interference correction □^^(£') in 
addition depends on physics at long distance scales |184l 11811 1182[ I200"| 1201] . The current best theoretical 
estimate from Ref. |201j is = 0.0713(8). An explicit energy- dependence is shown for □^^(ii^) in 
Eq. flllOp . in anticipation of a significant variation of this quantity for electron scattering in the GeV 
range. The other radiative corrections in Eq. flllOp are not expected to show such a dependence. 

Note that the absolute correction to the proton weak charge is related to the relative correction 
b^z in Eq. ^ by 



□,^(0) = Q%,b,z = Q'y^ 



TZeiM; 



W / _ (pv)\ 



:iiii 



Corrections from the interference of A^^^'' with the TPE amplitude under the replacement M.^ — > 
M.^ + M.^^ (namely, ^^(77) in Eq- (195!) ) vanish in the forward limit, and therefore do not affect the 
asymmetry. 

A precise determination of the proton weak charge can provide an important test of the Standard 
Model through verification of the predicted running of the weak mixing angle from the Z boson pole to 
low energies |201j . The Qweak experiment at Jefferson Lab [179] . which will complete data taking in 2012, 
was designed to measure to a higher level of precision than previously possible. In combination with 
constraints from atomic parity violation |202] . the Qweak measurement aims to either discover evidence 
for new physics beyond the Standard Model that leads to parity violation in electron scattering, or raise 
the limit on its mass scale to above 2 TeV, complementing direct searches at the LHC |20H 1203] . 

As discussed above, the correction in general has contributions from the vector electron-axial 
vector hadron coupling of the Z boson i^^z) ^'^'^ from the axial electron-vector hadron coupling of the 
Z (D^^^), 07^ = ^^z + '-'tZ- The vector hadron contribution D^^: vanishes in the limit of zero energy, 
but is finite at ii^ > 0. The axial hadron correction n^^, which is dominant at very low electron energies 
relevant to atomic parity-violation experiments, was estimated some time ago by Marciano and Sirlin 
|181[ 1182] in terms of a free quark model-inspired loop calculation. 

At forward angles, the correction can be computed from its imaginary part using standard 
forward dispersion relations analogous to Eq. f l53|) |187] . 

n.u.,(E^^\vf_jE'^-^^l^. ,112, 

The imaginary part of D-^^ depends on the PV ep — )■ eX cross section, which can be expressed in terms 
of interference electroweak structure functions. The integration over negative energies in Eq. f lll2p 
corresponds to the crossed 7Z box diagram, and can be related to the box diagram using the crossing 
symmetry properties of the 7Z amplitudes in Eq. (p2|) and (!93|) . For the vector hadron and axial hadron 
loop corrections one has, in analogy with Eq. (!54|) . 



9 771 /-CXD -I 

^en^^(i?) = —V \ dE'^^^^-^^XmuX^.iE'), (113) 



o /"OO 771/ 

^en^^(E) = -V\ dE'^j^^-^^XmU^,{E'). (114) 

'J 
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It is evident from Eq. flllSp that at zero energy the vector hadron correction vanishes, TZe □^^(0) = 0. 

From the optical theorem, the imaginary part of PV 7Z exchange amphtude can be written in terms 
of the cross section for all possible final hadronic states, 

2X™Kv'.A.r)^ -Me^(^) /^(^)rT^^;?-^". (115) 

where I = k — qi, Ei = ^J{k — qiY + m^, and Q\ = —qj. The interference hadronic tensor can be 
parametrized in terms of three interference electroweak structure functions, 

MW!:^ = -gt^'^F^^ + ^^F^^ -ie'^'^P^^^F^^. (116) 
^ P-Qi 2p ■ gi 

The structure functions F^2 ^^^^ analogous to the electromagnetic structure functions in Eq. ( IST!) and 
contribute to the vector hadron correction D^^, while the F^^ structure function contributes to the 
axial hadron correction n^^- They are in general functions of the exchanged boson virtuality and 
of the invariant mass W of the exchanged boson and proton, or alternatively of the Bjorken variable 
x = Ql/iW'-M^ + Ql). 

In analogy with Eq. ( 152|) . the imaginary parts of the vector hadron and axial hadron contributions 
to the 7Z box diagrams can be written as [Wl HHH M 12021 EDS] 



XmU^z{E) = Ve———: / dW 



ni /"Vl max rlt^'^ 



7Z I (Ql,max Ql) jp-yZ 

1 + Ql/Ml 

'3l,max rlD^ 



1 T^aTTTP). TTFTT^ 2 



2ME 1 



(117) 

F3^^ (118) 



^^^"^ ^^i2MEy Jm2'^" Jo l + Q\/Ml[W^-M^ + Ql 2_ 

where -Oe = (1 - 45^), with = sin^ 6 w{M^) = 0.23116(13) in the MS scheme [183]. Note that 
Eqs. flll7p and (11181) assume that the scale dependence of a and sin^ 9w is negligible (see, however, 
Ref. ^206j ). 

In evaluating the imaginary parts of it is convenient to split the integration into three regions: 
(i) elastic with W"^ = M^; (ii) resonance with (M + m^r)^ < < 4 GeV^; and (iii) deep inelastic 
scattering (DIS) region, with W"^ > 4 GeV^. The elastic contributions, including both the ■yZ and Z7 
orderings, can be written in terms of the elastic electroweak form factors defined in Sec. I7.H 

F^^^'''\W^QI) = 2M\GZGl5{W^ -M^), (119) 

Ff^^'\w\Ql) = ^ (gJG'/ + rGZGff) SiW - M'), (120) 

f;J^^'''\w^,QI) = -AM\GZGl5{W^ - M^), (121) 



where here r = Qj/4M^. Note that with dipole parametrizations of form factors the integrals flll3p . 
(11140 . (I117p and (I118P can be performed analytically, providing a useful cross-check of the numerical 
calculations. One can verify that the results in fact coincide exactly with the direct loop calculations 
[Bl 1185] of the 7Z corrections in Sec. 17.2. H which do not use dispersion relations, and in which the 
intermediate nucleon is off-shell. The result for TZe\^^2;{E) ^Iso agrees at = with the classic 
calculation of Marciano and Sirlin |181] . if the parameters are adjusted to correspond to those of 
Ref. |181] . To simplify notation, in the following we denote T^e by D^^^, since these are the 

quantities of interest for the observables. Numerical results for the elastic contribution, n'^^^\E), are 
given in Table [1] at several energies E. 

In the nucleon resonance region, while there is an abundance of electroproduction data, there are 
no direct measurements of the interference structure functions F^2 3- Fo'^ vector transitions to isospin 
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Table 1: Contributions to the real part of n^^[E) from different kinematical regions |2U6] . 
The quoted uncertainty on the total accounts for mo del- dependence as well as uncertainties 
in parameters. The energy of the Q„eak experiment is 1.165 GeV. 







Hi V j 






0.0 


1.165 


3.0 


elastic 


0.00064 


0.00005 


0.00001 


resonance 


0.00023 


0.00012 


0.00002 


DIS > 1 GeV^) 
DIS < 1 GeV^) 


0.00327 


0.00328 


0.00330 


0.00024 


0.00025 


0.00027 


Total 


0.0044(4) 


0.0037(4) 


0.0036(4) 



/ = 3/2 states, such as the A(1232) resonance, conservation of the vector current and isospin symmetry 
require the weak isovector transition form factors to be equal to the electromagnetic ones multiplied 
by (1 + Q^). For isospin J = 1/2 resonances, which contain contributions from isovector and isoscalar 
currents, using SU(6) quark model wave functions one can verify that for the most prominent I = 
1/2 states the magnitudes of the Z'-boson transition couplings are equal to the respective photon 
couplings to within a few percent |189[ 1204] . Estimates of the resonance region contribution to the 
box diagram, n]^^'^^\ have been made recently |187[ I189[ I204[ 1205] by parametrizing resonance region 
structure function data from SLAG, Jefferson Lab and elsewhere in terms of resonant and nonresonant 
background components |204[l207j . The model dependence of relating the vector jZ resonance structure 
functions to the electromagnetic structure functions has been studied in Ref. [189] . 

For the axial- vector resonance contributions, n^^^^\ one can use parametrizations of the transition 
form factors in neutrino scattering from Lalakulich et al. \195\ 1196] , with modified isospin factors 
appropriate to 7Z. These form factors have been fit to the Jefferson Lab pion electroproduction data 
(vector part) and pion production data in and u scattering at ANL, BNL and Serpukhov (axial- 
vector part) up to Ql = 3.5 GeV^. At larger the resonance contributions are suppressed by the Ql 
dependence of the transition form factors, which is stronger for the dominant A(1232) resonance than 
for the higher- mass resonances |196] . The resulting resonance contribution, n^^'^^\o), is smaller than 
the elastic term at E = 0, but decreases less rapidly with increasing energy, as shown in Tabled] 

The contributions from high W can be computed by analyzing the high-Q^ (Ql > Qq) and low-Qi 
{Ql < Ql) regions separately, with Ql taken to be ~ 1 GeV^. Structure functions in the high-VT, 
high-Qi region of deep-inelastic scattering can be related to leading twist parton distribution functions 
(PDFs), and at leading order in as are given by 

^.z(Dis) ^ ,J22e,gU^ + Q) = 2xFf^^''\ (122) 

^^,z(Dis) ^ a;^2e,^l(g-g), (123) 
q 

where q and q are the quark and antiquark PDFs, and gy = 1/2 — (4/3) sin^ 9w and Qy = —1/2 + 
(2/3) sin^ Ow are the weak vector charges for u and d quarks, respectively. Note that in the limit 
2gy — Cq the interference structure functions F^^^'^^^^ — > F2^^^^\ where 

F^^^) = xJ2el{q + q). (124) 
In the case of a free quark target one has F^l^ = (5/3) x 6{1 — x), which gives a contribution to the axial 
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hadron 7Z correction 

^-^^ = Ve - (ln^ + -). (125) 

This reproduces the perturbative result in the free quark model of Ref. [181] at = 0, where M is 
interpreted hadronic mass scale. 

As observed in Ref. |187] . for instance, for three quark flavors the sum over electroweak charges 
J2q'^^q9v ~ §(-'- + Qw) ~ Tlq^^ so that in the high-W^ region dominated by sea quarks, the quark 
distributions are approximately flavor independent, in which case F^f ~ -F72- addition, while 
the Callan-Gross relation between the Fi and F2 structure functions is assumed in Eq. (11221) at high 
Q\, at finite Q\ the contribution from Fi is obtained from F2 and the ratio Rlt = ctl/c^t = (1 + 
4M^x^/Q^)-F2/(2xFi) — 1 of the longitudinal and transverse cross sections |208] . 

One can simplify the expressions for nl^'z in the DIS region by firstly interchanging the order of the 
integrations in Eqs. (I113p . (I114p and (I117p . (11180 to perform the integrals over the energy analytically 
[205(1206] . The resulting integrands can be expanded at high Q\ and low E in powers of /Q\, yielding 
a series whose coefficients are moments of structure functions [206j . 



= -2ME r 

^ Jo 



Q\{1 + Ql/Ml) 



^^(4E2_5g2)^(4)^__ 



5Q1 



^ 3a r dQ^^ 



Here the moments of the structure functions are defined as 



Min^(5E^-3Q?)Mf + 



(126) 
(127) 



M^-\Ql) ^ [\xx-'j^f{x,Ql), ^ = 1,2,3, 
^0 



128) 



where J-^^ = jxF^^, Fg*^^, xFg Note that the upper limit Xmax on the x-integrals in Eqs. (I126p and 
(I127p has been approximated by 1; the resulting error was found to be less than 10"'' for Qf > 1 GeV^ 
|. The large-x contributions to M^"'\qI) become more important for large n; however, the higher 



moments are suppressed by increasing powers of l/Qf- In practice, the integrals in Eqs. (11261) and 
(I127p are dominated by the lowest moments, with the l/Qf corrections being relatively small in DIS 
kinematics. For the axial- vector hadron part, the lowest moment, mI^^\qI), is the 7Z analog of the 
Gross-Llewellyn Smith (GLS) sum rule [209j for uN DIS, which at leading order counts the number 
of valence quarks in the nucleon. The corresponding quantity for jZ is 26,(7^ = 5/3, so that at 
next-to-leading order in the MS scheme 



Ml'\Qi) = (129) 



The lowest moment contribution to Eq. (I127p is therefore 



3« r ,^2 M^^\Ql 
'2^JqI ^'Ql{i + Qi/Mh 



□T'HO) - / dQl—p-^lL^, (130) 



which is identical to the Marciano-Sirlin result [181] for the high energy contribution to the box diagram. 

For Ql < Qq a partonic description of the structure functions is not valid. In particular, since the 
integrals over Ql in Eqs. (I117P and (I118P extend down to Qi = 0, and the upper limit on the x-integral. 
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^^max, is also limited by Q^, one requires the behavior of the structure functions in the limit of both low 
X and low Qi- In the case of the vector structure function, conservation of the two vector currents 
requires ~ Ql as Ql — )■ 0. In computing the high-iy contributions to several authors have 
used Regge inspired parametrizations |123[ I210[ 1211] of the electromagnetic structure functions, and 
approximating ^ F2 at small x. 

By contrast, F^^ depends on both vector and axial- vector currents, so that no analogous current 
conservation constraint exists. In the absence of data on F^^ in the low-x, low-Qi region, Blunden et 
al. considered models for its possible x and Q\ dependence such that F^'^ at Xmax should not diverge 
as Q\ — 7- 0, and should match the partonic structure function at Q\ = Q^. This model-dependence is 
reflected in the quoted uncertainty in Table [1] 



1.2 r 



o 





E (GeV) 

Figure 39: 7Z box corrections to for the (left) vector hadron contribution show- 
ing the resonant (dashed) and nonresonant (dotted) components, and the sum (solid, and 
shaded); (right) axial hadron contribution (labeled "A") together with the sum of axial 
and vector hadron corrections ("V+A"), and the F = result from Refs. (TEH [200] ("MS", 
extended to finite E for comparison). The vertical lines ai E = 1.165 GeV indicate the 
energy of the Qweak experiment. Figures taken from Refs. )204] 1206] . 



The total vector and axial hadron corrections (E) are shown in Fig. [39] as a function of the 
incident electron energy E. For the vector hadron correction D^^, Fig. [39] (left), most of the strength 
(~ 80%) comes from relatively low energies, below 4 GeV, where the Qi range extends to ~ 6 GeV^, 
and to ~ 3 GeV. This coincides precisely with the region that a wealth of very accurate electro- 
production data exists from Jefferson Lab )146] 1212] 1213] 1214] . The nonresonant contribution to 
is small at low energies, rising linearly with E in this region. The resonant part increases steeply to 
a maximum at F ~ 1 GeV, before falling off like 1/E |187l 1204] . Sibirtsev et al. find the resonant 
and nonresonant contributions to to be 0.0026 and 0.0021, respectively, at the energy relevant for 
the Qweak experiment, E = 1.165 GeV (although one should note that this separation is somewhat 
arbitrary, as only the total cross section is physically meaningful). The error band, which grows with 
energy, is obtained from the uncertainty in the fit parameters using a variational method |204] . 

Using different structure function inputs, Rislow and Carlson |205] find □^2(1-165 GeV) = 0.0057 ± 
0.0009, which is slightly higher than but consistent with the value from Sibirtsev et al. )204] . Gorchtein 
et al. H89j, on the other hand, find n;^^ (1.165 GeV) = 0.0054 ± 0.0020, which is again consistent with 
the other estimates but has a larger uncertainty due to the larger range of input structure functions 
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considered there. 

The axial hadron correction in Fig. (right) is dominated by the DIS contribution, which has 
neghgible E dependence. On the other hand, the resonance and low-Q^ DIS contributions dominate 
the uncertainties. The total axial hadron correction 0^^(^E) is 0.0044(4) at E = 0, and 0.0037(4) at 
E = 1.165 GeV. This should be compared to the value 0.0052(5) used in Ref. |200] . which is assumed 
to be energy independent. Combined with the correction to D^^, this shifts the theoretical estimate for 
from 0.0713(8) to 0.0705(8), with a total additional energy dependent correction of 0.0040]^q'qqq4 at 
E = 1.165 GeV. 
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Figure 40: Dependence of the vector hadron correction on the momentum transfer 
squared t at the Qweak energy E = 1.165 GeV. The outer curves indicate the uncertainty. 
Figure taken from Ref. jl89] . 



The t dependence of the corrections, from the forward limit — t = = to the Qweak value of 
—t = 0.03 GeV^, was considered by Gorchtein et al. |187l 1189] using the phenomenological ansatz 

n^.iE, t) = □,^(E, 0) "^P^~f^f^'/^\ (131) 

where F^^{t) is the Dirac proton form factor, and B = (7 ± 1) GeV~^ is extracted from analysis of 
data forward Compton scattering data on ^He nuclei |215] . Although the estimate is considered by the 
authors to be exploratory, it suggests that the dispersion correction decreases by only about 2% from 
its t = value [189J. 

The corrections D^^^ are important for the interpretation of the Qweak experiment, given its projected 
uncertainty of ±0.003 |179j . It is also critical to the physical interpretation of the experiment which is 
expected to constrain possible sources of parity violation from beyond the Standard Model at a mass 
scale of > 2 TeV |203] . The uncertainties in the corrections can be reduced with future parity- violating 
structure function measurements at low Q^, such as those planned at Jefferson Lab [216^ 1217] . The 
high precision determination of Q^r would then allow more robust extraction of signals for new physics 
beyond the Standard Model. 



8 Conclusions and outlook 

The renewed interest in the role of two-photon exchange in electron-hadron scattering sprouted from 
the challenge to resolve a major discrepancy between new polarization transfer measurements of the 
proton's electric to magnetic form factor ratio and earlier results using Rosenbluth separations of elastic 
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cross sections, which has been confirmed now by several independent experiments. In the intervening 
decade since the first polarization transfer experiments were performed at Jefferson Lab, the study of 
TPE has flourished into an area rich in phenomenology, with many developments in experiment and 
theory that have significantly altered how we think about the basic process of electromagnetic scattering 
from hadrons. The most likely candidate that has emerged to explain the form factor discrepancy is 
two-photon exchange, in particular the hadron structure dependent parts of the 77 box and crossed-box 
diagrams neglected in the standard treatments of radiative corrections. 

In this article we have reviewed the developments in the study of TPE over the past decade, focusing 
on the hadronic framework for radiative corrections applicable for most measurements extending into 
the ~ few GeV^ range. This period has witnessed major advances in theoretical calculations of TPE 
amplitudes, which we have outlined in this review, together with their implications for observables. 
We have surveyed the effects of TPE corrections on various elastic ep scattering observables, including 
Rosenbluth separations of cross sections, polarization measurements, and e^p/e~p ratios, as well as the 
Born-forbidden normal beam and target asymmetries. Two-photon effects also find their way into other 
observables, such as neutron form factors, resonance electroproduction, the pion form factor, and elastic 
electron-nucleus (deuteron and ^He) cross sections, which we have summarized here. The theoretical 
advances in TPE computations have furthermore allowed for the first time the a priori inclusion of TPE 
corrections in a global fit of form factor data, which has recently been performed. 

Extending the discussion into the weak sector, we reviewed the calculation of corrections to parity- 
violating elastic scattering cross sections and asymmetries from 7Z interference diagrams. These turn 
out to have relatively minor impact on the extraction of the strange nucleon electromagnetic form 
factors, but play an important role in measurements of the weak charge of the proton at forward angles, 
which gives direct access to the weak mixing angle, sin^ 6w- 

At present, the strongest experimental constraints on the size and kinematic dependence of TPE 
corrections come from the discrepancy between form factors extracted from elastic cross section and 
polarization measurements, assuming one attributes the discrepancy to TPE. Measurements of beam 
normal single spin asymmetries also provide important tests of TPE calculations, but are sensitive to 
the imaginary parts of TPE amplitudes which are not directly probed in form factor measurements. 
The most direct experimental information comes from comparisons of positron and electron scattering 
from the proton at large e, and the limit on nonlinear contributions in e to polarization observables. 
Unfortunately, at present these provide only limited evidence for TPE contributions, and serve mainly 
to constrain the magnitude of the effects. 

As we have shown in this review, modern calculations of TPE corrections can resolve most of the 
observed form factor discrepancy and are consistent with all other experimental constraints; however, 
it will be important to test their validity over as large a range of observables and reactions as possible. 
Including the TPE corrections and applying an estimate of their uncertainties at high yields form 
factors in which TPE effects are not the dominant source of uncertainty |129] . 

In future an important goal will be to obtain direct experimental verification that the form factor 
discrepancy is indeed explained by TPE effects. This will involve an experimental program consisting 
of several efforts, including measurement of the e^p/e~p cross section ratio and detailed examination of 
the e dependence of polarization observables, particularly since in some cases theoretical predictions for 
these vary considerably. An important opportunity may be presented by the development of positron 
beams at a future 12 GeV upgraded Jefferson Lab [218j . In addition, measurement of various beam and 
target normal asymmetries will provide direct evidence for Born-forbidden effects and constrain input 
for dispersion relation analyses. 

Further theoretical work will seek better control of higher-mass intermediate state contributions to 
the box diagrams, especially at higher values. Along this direction, it will be important to go beyond 
the resonance approximation to include the nonresonant background in the intermediate state spectrum. 
This naturally becomes difficult to do within any specific model, and a more effective approach may 
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be to utilize dispersion relation to obtain TPE amplitudes from the imaginary parts of the Compton 
scattering data. This is particularly appealing at forward angles, where the amplitudes can be related 
to inclusive structure functions. 

This approach is proving to be especially relevant for estimating 7Z interference corrections to the 
proton weak charge at low Q^, currently being measured in the Qweak experiment at Jefferson Lab. 
Extending this program into the non-forward region will be an important challenge, but one which 
will allow for a reliable global analysis of the parity-violating elastic scattering data to extract strange 
nucleon form factors correcting for the effects of two-boson exchange. 

Finally, we should note that TPE corrections are indeed small in general, and for a large range of 
experiments are at the 0{a) (or ~ 1%) level that has historically been assumed. It is only certain 
cases, typically at backward angles for two-photon exchange or forward angles for 7Z corrections, 
where they tend to be larger, at the ~ few % level; is a very special case where even a few percent 
effect in the cross section has a large impact on the Rosenbluth Ge/Gm extraction. Two-photon 
exchange corrections therefore do not necessarily require revising methodologies for electron scattering 
in general, so much as reminding us of the need to take greater care in estimating their possible impact 
on observables that may be particularly susceptible to their effects. 
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